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ABSTRACT 



Aims. We aim to study the properties of massive stars at low metallicity, with an emphasis on their evolution, rotation, and surface 
abundances. We focus on O-type dwarfs in the Small Magellanic Cloud. These stars are expected to have weak winds that do not 
Ci remove significant amounts of their initial angular momentum. 

Methods. We analyzed the UV and optical spectra of twenty-three objects using the NLTE stellar atmosphere code CMFGEN and 
I derived photospheric and wind properties. 

Results. The observed binary fraction of the sample is ss 26%, which is consistent with more systematic studies, if one considers 
that the actual binary fraction is potentially larger owing to low-luminosity companions and that the sample was biased because it 
excluded obvious spectroscopic binaries. The location of the fastest rotators in the Hertzsprung-Russell (H-R) diagram built with 
fast-rotating evolutionary models and isochrones indicates that these could be several Myr old. The offset in the position of these 
fast rotators compared with the other stars confirms the predictions of evolutionary models that fast-rotating stars tend to evolve more 
vertically in the H-R diagram. Only one star of luminosity class Vz, expected to best characterize extreme youth, is located on the 
ZAMS, the other two stars are more evolved. We found that the distribution of O and B stars in the e(N) - v sin i diagram is the same, 
which suggests that the mechanisms responsible for the chemical enrichment of slowly rotating massive stars depends only weakly on 
the star's mass. We furthermore confirm that the group of slowly rotating N-rich stars is not reproduced by the evolutionary tracks. 
Even for more massive stars and faster rotators, our results call for stronger mixing in the models to explain the range of observed N 
abundances. All stars have an N/C ratio as a function of stellar luminosity that match the predictions of the stellar evolution models 
well. More massive stars have a higher N/C ratio than the less massive stars. Faster rotators show on average a higher N/C ratio than 
slower rotators, again consistent with the expected trend of stronger mixing as rotation increases. When comparing the N/O versus 
N/C ratios with those of stellar evolution models, the same global qualitative agreement is reached. The only discrepant behavior 
is observed for the youngest two stars of the sample, which both show very strong signs of mixing, which is unexpected for their 
r- H evolutionary status. 

t> Key words. Stars: early-type - Stars: fundamental parameters - Stars: rotation - Stars: abundances - Galaxies: Magellanic Clouds 
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1 . Introduction eter of stellar evolution models because O-type stars loose a siz- 
able fraction of their total mass during their lifetime. Progress 

During the past three decades, much effort has been invested has been achi eved in the parameterization of mass-loss rates 

in understanding the effect of essential processes such as stel- on basic physical parameters, such as luminosity and metallic- 

lar winds (and the resulting mass loss) and stellar rotation on the ity> from the the0 retical and the observational sides (Vink et al. 

evolution of massive stars. Mass-loss rates are a crucial param- 2 001; Krticka 2006; Mokiem et al. 2007; Muijres et al. 2012). 

Although the general paradigm is successful in explaining the 

Send offprint requests to: J.-C. Bouret global behavior of stellar winds, several important questions 

* Based on observations made with the NASA-ESA Hubble Space rema in open, including the presence of structures in the wind 

Telescope (program GO 1 1625), obtained at STScI, which is operated (clumping) and the so-called weak winds in late-O dwarfs (see 

byAURA Inc. under NASA contract NAS 5-26555 Puis etal 2008, for a review). 



Based on observations collected at the European Southern Obser- 
vatory Very Large Telescope, program 079.D-0073 
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Likewise, major progress in stellar rotation theory has been 
achieved in the past decade. The effects of rotation include an 
increase of the mass-loss rate, a change in evolutionary tracks in 
the H-R diagram, a lowering of the effective gravity, an exten- 
sion of the main-sequence phase, and the mixing of CNO-cycle 
processed material up to the stellar surface (Maeder & Meynet 
2000; Heger & Langer 2000). Qualitatively, stellar evolution 
models that include rotation predict a surface enrichment of he- 
lium and nitrogen with an associated carbon and oxygen deple- 
tion during the main-sequence evolution. Quantitatively, stellar 
evolution models with rotation predict that fast-rotating stars un- 
dergo more mixing, and hence a greater nitrogen (and helium) 
surface enrichment should be observed. Because of the sensi- 
tivity of evolutionary tracks in the H-R diagram with rotation, 
the location of massive main-sequence stars in the H-R diagram 
does not provide a unique determination of age and mass. Ro- 
tation thus introduces some degeneracy, but surface abundance 
patterns could help resolve this degeneracy (Brott et al. 2011a; 
Kohler et al. 2012). 

Hot-star winds are driven by momentum transfer from the 
radiation field to the gas via photon absorption in spectral lines 
of metal elements. A dependence of the mass-loss rates on the 
stellar metallicity is thus predicted (M oc z°- 69±0lQ - Vink et 
al. 2001) and observed (M oc z a83±016 - Mokiem et al. 2007). 
In addition to mass, the stellar winds take angular momentum 
away, and therefore the removal of angular momentum during 
the star's evolution is significantly reduced at lower metallici- 
ties, increasing the effects of rotation (Maeder & Meynet 2001). 
The interplay between stellar rotation and the stellar wind, which 
blurs the stars' initial rotational velocity properties within only a 
few million years in the case of Galactic massive stars, is thus ex- 
pected to be reduced with decreasing metal content. The initial 
conditions of rotational velocity remain better preserved during 
the main-sequence life of O-type stars at low metallicity. Mas- 
sive stars in the Magellanic Clouds (MCs) are therefore the ap- 
propriate targets to study the effects of rotation on the evolution 
of massive stars. 

Several studies have tackled this topic, though from different 
perspectives. From the analysis of four extreme O-type super- 
giants, Crowther et al. (2002) found that the CNO abundances 
differ greatly from those inferred from nebular and stellar stud- 
ies in the MCs. Nitrogen was strongly enhanced, with carbon 
(and oxygen) moderately depleted, which is suggestive of mix- 
ing of unprocessed and CNO-processed material at their sur- 
faces. These results were supported by Hillier et al. (2003), 
who focused on two giant O stars in the Small Magellanic Cloud 
(SMC) that have very similar effective temperatures and lumi- 
nosities, but whose wind spectral signatures differ markedly. Ni- 
trogen is substantially enhanced at the surface of one star, while 
the N/C abundance ratio in the other is below solar which agrees 
with SMC nebular abundances. A consistent picture of the evo- 
lution of these two stars may be obtained if they both had a simi- 
lar initial mass, M « 40 M Q , one star being a slow rotator and the 
other a fast rotator. Bouret et al. (2003) analyzed six O dwarfs in 
the SMC giant H II region NGC 346 and found that a majority of 
main-sequence O stars show enhanced nitrogen at their surface, 
indicating either fast stellar rotation and/or very efficient mixing 
processes. 

The Hillier et al. (2003) and Bouret et al. (2003) results were 
extended to the complete sample of 18 O stars in the SMC (HST 
GO Program 7437, PI: D. Lennon) in Heap et al. (2006). The 
program SMC O-type stars were chosen for their narrow-lined 
spectra to avoid problems with line blending in the analysis. 
Heap et al. (2006) found that the surface composition of a high 



percentage of the stars (about 80% of the sample) show surface 
nitrogen enrichment during their lifetime on the main sequence. 
Furthermore, there is no enrichment in helium and no depletion 
in carbon or oxygen, although the latter is poorly constrained, 
because we lack reliable abundance diagnostics. N/C enrich- 
ment factors of 10 - 40 were found, which is significantly higher 
than the theoretical predictions of factors of 5 - 10 (Maeder & 
Meynet 2001). In the few cases where nitrogen surface abun- 
dance has not yet been affected by mixing, Heap et al. (2006) 
also obtained CNO abundances consistent with SMC nebular 
abundances. Most stars in their sample show the so-called mass 
discrepancy between spectroscopic and evolutionary mass esti- 
mates. Along with nitrogen enrichment, Heap et al. (2006) in- 
terpreted this discrepancy as a result of rapid rotation that lowers 
the measured effective gravities. Therefore, most stars in their 
sample should be fast rotators with v rot > 200 kms -1 , the low 
measured u sin i thus implying that most stars are viewed almost 
pole-on (note again that their sample was biased toward stars 
with sharp spectral lines, hence low usin i). 

More recently, Hunter et al. (2007, 2008, 2009) investigated 
the link between surface abundance patterns and rotation on a 
larger, unbiased (against rotation) sample of Galactic and MC 
B-type stars observed for the VLT-FLAMES survey (Evans et al. 
2005, 2006). They concluded that up to 40% of their sample had 
nitrogen abundances that could not be explained by the theory of 
rotational mixing as implemented in stellar evolution models. In 
particular, they found populations of non-enriched fast rotators, 
highly enriched slow rotators, and B-supergiants that are highly 
enriched compared with normal core-hydrogen-burning objects. 

The results above have been used to calibrate the amount 
of convective overshooting and the efficiency of rotationally in- 
duced mixing in stellar evolution models (Brott et al. 2011a). 
More observational constraints were obtained using single-star 
population-synthesis and an extended grid of evolution mod- 
els, to simulate population of stars with masses, ages, and rota- 
tional velocity distributions consistent with those from the VLT- 
FLAMES sample (Brott et al. 201 lb). Although the observed 
fraction of stars without significant nitrogen enrichment could be 
reproduced, the simulations also predicted too many rapid rota- 
tors with enhanced nitrogen. Other problems included the failure 
to account for populations of slowly rotating, nitrogen-enriched 
objects on one hand, and of rapidly rotating non-enriched objects 
on the other hand. Brott et al. (201 lb) concluded that in addition 
to rotational mixing, other mechanisms should be accounted for, 
including the effects of binarity and magnetic fields on stellar 
evolution. 

It is therefore fair to state that the relative influence and in- 
terplay of mass loss, metallicity, age, rotation, and binarity on 
the global evolution of massive stars are not yet fully understood 
(Maeder et al. 2009). To answer these questions, we must ac- 
curately place stars on the theoretical H-R diagram and examine 
their other properties, such as rotation rates, masses, mass-loss 
rates, and surface CNO abundances. Many of these quantities 
can be accurately determined from spectral diagnostics in the 
UV range. This spectral range contains all major UV resonance 
wind lines and numerous other metal lines, providing the essen- 
tial diagnostics of the photospheric and wind properties of the 
stars (see e.g. Hillier et al. 2003; Bouret et al. 2003). 

We present in this paper the first results of an extensive anal- 
ysis of the UV spectra of 23 O-type dwarf stars in the SMC, 
with a special emphasis on the evolution and chemical proper- 
ties. In Sect. 2, we introduce the HST observing program, and 
the additional optical data used in this work. Sect. 3 presents the 
modeling tool and strategy, while general results are presented in 
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Sect. 4. We discuss the properties of the sample stars in Sect. 5 
and summarize in Sect. 6. Atomic data used in the modeling 
are disccused in Appendix A. Properties of individual stars are 
presented in Appendix B. 



2. Target selection and observations 

2.1. General presentation of the FUV program 

We aimed at mapping all spectral types and luminosity classes in 
the O-type range. For this aim, we secured high-resolution far- 
ultraviolet (FUV) spectra of 22 O-type stars in the SMC with the 
Cosmic Origin Spectrograph (COS) onboard the Hubble Space 
Telescope (HST Program GO 11625). Combined with a sec- 
ond set of FUV spectra obtained earlier with HST/STIS (Pro- 
gram GO 7437, PI: D. Lennon), the whole sample contains 39 
O stars. The sample was selected mainly from the UBVR CCD 
photometric survey of the Magellanic Clouds by Massey (2002). 
This survey is photometrically complete to V « 15.7, although 
some stars could have been missed in crowded regions. Ta- 
ble 6 in Massey (2002) provides cross-identifications with previ- 
ous prism-objective surveys. We also added four O dwarfs that 
are members of the NGC 346 cluster and are included in the 
ESO/VLT-FLAMES survey by Evans et al. (2006). 

In a few instances, we included two or three stars of the same 
spectral type and luminosity class to directly compare stars that 
were presumably similar, but had different apparent rotational 
velocities (as indicated by spectral qualifiers, for instance "n"). 
In addition to investigating the effects of rotation, this approach 
provides a direct way to distinguish between general trends and 
individual anomalies. Rotation velocity, however, was not a se- 
lection criterion per se because of the lack of prior information - 
but for the spectral qualifier "n' " - for many selected stars. While 
the distribution of v sin i for OB stars in the SMC peaks around 
150-180 km s (Mokiem et al. 2006), only one star has v sin i > 
120 km s _1 (AzV 80) in the STIS sample (Heap et al. 2006). This 
sample therefore shows a preponderance of stars with narrow- 
lined spectra, which may suggest a bias toward slower rotators 
(Heap et al. argued, however, that the sample was biased toward 
stars seen mostly pole-on). Our selection of the COS sample 
aimed at offsetting this potential bias. We stress that the profiles 
of UV wind lines retain the information on wind properties even 
for fast-rotating stars with broadened line profiles, because of the 
high terminal velocities of O-star winds. 

Finally, color excesses and bolometric corrections were de- 
rived from the Massey photometry and from theoretical colors 
(Lanz & Hubeny 2003). From this list, we selected stars with 
the lowest extinction to maximize the expected FUV flux, and in 
turn the expected signal-to-noise ratio (S/N) in a given UV ex- 
posure. For each selected target, we checked also that there was 
no other bright nearby star that might be in the COS aperture. 

In this paper, we focus on 23 stars that are classified as 
dwarfs, because they are presumably closer to the ZAMS and 
less evolved than the giants and supergiants. One possible ex- 
ception is MPG 355, classified as ON2III (f*) by Walborn et al. 
(2004) but also as 02V by Massey et al. (2009). Taken at face 
value, this indicates that a non-ambiguous identification of lumi- 
nosity classes from spectral morphology becomes difficult for 
early-O-type stars. Problems with the luminosity class crite- 
rion (based on the morphology of He n /14686) may also be- 
come more severe at low-Z environments such as the SMC (see 
the very interesting discussion about the specific case of early- 
O-type stars in low-Z environments in Rivero Gonzalez et al. 
2012). Therefore, some early-O-type stars whose physical prop- 



erties (visual luminosities and surface gravities) suggest they are 
giants, might actually be dwarfs according to their (spectroscop- 
ically determined) luminosity classes. Overall, dwarf and giant 
stars occupy nearly the same location in the H-R diagram at the 
hot end of the main sequence (which is known to be the case for 
MPG 355 from Bouret et al. 2003). Accordingly, we elected to 
include MPG 355 is our study which will thus sample the full 
range from the earliest to the latest O stars. 

O-type stars spend most of their lifetime in the dwarf phase, 
while later stages are much shorter. We intend to investigate 
the effect of rotation on the surface abundances in this phase. 
Their stellar winds are also weaker. An analysis of the more 
luminous objects (luminosity classes III and I) will be presented 
in a separate paper, with a stronger emphasis on mass loss. 

2.2. UV observations and optical data 

We observed each target in a sequence of four science exposures 
with HST/COS, using gratings G130M and G160M with central 
wavelengths (A1291, ,11327) and 011577, /H623) for the adopted 
settings respectively. This sequence provides coverage of the 
FUV spectrum without a gap between 1132 and 1798 A. Over- 
lapping segments were co-added to improve the final signal-to- 
noise ratio. Thanks to the high efficiency of COS in the FUV, 
S/N « 20 - 30 per resolution element were typically achieved in 
a single HST orbit for each star. For SMC targets, the visibility is 
60 min per orbit. Accounting for various observation overheads, 
only 29 min were left for the science exposures. We therefore 
allocated 300 sec for the G130M exposures and 570 sec for the 
G160M exposures in order to achieve the best possible S/N over 
the whole spectral range. 

Some stars were observed previously with HST/STIS 
through the 0".2 - 0".2 aperture, using the far-UV MAMA de- 
tector in the E140M mode. A spectral interval from 1150 A to 
1700 A is covered in a single exposure at an effective spectral re- 
solving power of R = 46,000. With one exception, two or more 
observations per target were obtained. The exposure times were 
set according to the apparent magnitudes of the stars to equalize 
the S/Ns of the spectrograms. The total exposure times ranged 
from 45 minutes for MPG 355 (one observation) to 245 minutes 
for MPG 012 (five observations). At 1300 A (the wavelength 
with the highest sensitivity), the S/N per binned data point ranges 
from 55 to 110; at 1600 A it ranges from 25 to 50. A more de- 
tailed description of the STIS observations and reduction can be 
found in Walborn & Howarth (2000) and Heap et al. (2006). 

We supplemented the FUV spectra with optical spectra for 
most stars. First, we obtained spectra with the UVES spectro- 
graph on the VLT (Dekker et al. 2000) under program 079. D- 
0073(A). The spectra have a resolving power of about 45,000, 
with five pixels per resolution element. The data were reduced 
using the UVES context within MIDAS (see Ballester et al. 
2000). The reduction process includes flat-fielding, bias and sky- 
substraction, and a relative wavelength calibration. Cosmic-ray 
removal was performed with an optimal extraction method for 
each spectrum. An additional VLT/UVES spectrum was ob- 
tained for AzV 388 in program 074.D-0164, and was kindly 
provided by P. Crowther. The optical spectrum of AzV 177, 
obtained at CTIO and presented in Massey et al. (2009), was 
kindly provided by P. Massey. For those stars with HST/STIS 
spectra, optical spectra were available from observations with 
AAT/UCLES or ESO/CASPEC. They were extensively pre- 
sented in Walborn et al. (2000) and Heap et al. (2006), to which 
we refer for more details. Optical spectroscopy was also supple- 



Article number, page 3 of 43 



mented from the ESO VLT/FLAMES survey of massive stars in 
the SMC (Evans et al. 2006). See Evans et al. (2005, 2006) for 
a general description of the properties of the FLAMES data and 
their reduction. For a few more objects (see Table 1), we used 
lower-resolution spectra obtained for the AAT/2dF program in 
the SMC in Evans et al. (2004). These authors provide details 
about the data and their reduction. In total, we secured optical 
spectra for 18 stars out of the 23 stars of our sample. 

The optical coverage is especially important in determining 
the stellar surface gravities, because this can only be poorly con- 
strained from the UV spectrum 1 . However, optical spectroscopy 
presents us with some particular difficulties in determining the 
stellar parameters, in particular because of nebular emission fea- 
tures in optical spectra (Heap et al. 2006; Mokiem et al. 2006). 
This nebular masking, for the Balmer and He i lines that are 
used to constrain T e g, is especially critical for some stars in 
NGC 346. Although the nebular features are clearly resolved 
in VLT/FLAMES spectra, an efficient and accurate substraction 
from photospheric line cores remains problematic (Evans et al. 
2005). It is difficult to assess the actual level of filling that is 
left after removal of the nebular emission, and the He i lines may 
seem weaker than they actually are, which would bias the T e ff 
estimates toward higher values. This problem is also relevant for 
field stars observed at lower resolution (2dF spectra). We dis- 
cuss our methodology of determining effective temperatures and 
surface gravities in more details in §3.2. 

3. Spectroscopic analysis 

3.1. Model atmosphere 

We performed the analysis using the model atmosphere code 
CMFGEN (Hillier & Miller 1998; Hillier et al. 2003). This code 
computes NLTE line-blanketed model atmospheres, solving the 
radiative transfer and statistical equilibrium equations in the co- 
moving frame of the fluid in a spherically-symmetric outflow. To 
facilitate the extensive inclusion of line blanketing, a formalism 
of super-levels was adopted. This approach allowed us to incor- 
porate many energy levels from ions of many different species 
in the model atmosphere calculations. More specifically, we in- 
cluded in this work ions of H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, 
and Fe, which represent up to « 14,000 individual levels. After 
convergence, a formal solution of the radiative transfer equation 
was computed in the observer's frame (Busche & Hillier 2005), 
thus providing the synthetic spectrum for a comparison to obser- 
vations. 

The photospheric density structure was initially described 
with an hydrostatic density structure from a Tlusty model at- 
mosphere (Hubeny & Lanz 1995a; Lanz & Hubeny 2003) in the 
quasi-static deep layers, while we used a standard/?- velocity law 
in the wind. The /^-velocity law is connected to the hydrostatic 
density structure just above the sonic point (at approximately 15 
kms" 1 ). The mass-loss rate, density, and velocity are related via 
the continuity equation. 

In more advanced models we iterated on the density structure 
so that the hydrostatic equation, 



dP/dr = -g gmv + g lad , 
where 

= GM/r 2 



9: 



grav 



(1) 

(2) 



1 It is well known that UV wind features for O stars correlate with lu- 
minosity class (e.g. Walborn & Panek 1984; Walborn & Nichols-Bohlin 
1987). Because we did not derive M and v^ from first principals, we 
cannot use this effect to further constrain the photospheric parameters 



and the radiative acceleration g m a 

4tt r 
gmi = — P I XvH v dv (3) 

was satisfied. We solved for the density structure by integrat- 
ing the hydrostatic equation using the Runge-Kutta method. To 
facilitate the integration, the Rosseland LTE opacity was com- 
puted as a function of temperature and electron number density 
for the model abundances prior to the CMFGEN run. In a typ- 
ical CMFGEN model, a hydrostatic iteration is performed after 
the first model iteration and then every n (with n ~ 8) iterations. 
Generally three to four hydrostatic iterations are needed to ob- 
tain convergence to better than 5% everywhere. 

Model atmospheres were calculated assuming that the winds 
were clumped, for all stars showing obvious spectral signatures 
of a wind (i.e., P Cygni line profiles). We adopted a paramet- 
ric treatment of optically thin clumps with a void interclump 
medium that is described by two parameters, a volume filling 
factor / and the velocity where clumping starts v c \. This sim- 
ple approach has some limitation such as a monotonic descrip- 
tion, while more elaborate theoretical simulations by Runacres 
& Owocki (2002); Puis et al. (2006) instead suggest that clump- 
ing is non-mono tonic. The radial distribution of the clumping 
factor was constrained in Puis et al. (2006) and Najarro et al. 
(2011), from simultaneous modeling of Ha, infrared, millime- 
ter, and radio observations. Nevertheless, the actual variation 
is still uncertain, and we therefore adopted this reasonable sim- 
ple approach. When no clear wind signatures are observed, we 
restricted ourselves to calculating homogeneous wind models. 

Since Mhom - -^ciump/ v/' tne mass -l° ss rates that we derived 
must be considered potentially as upper limits (and these mass- 
loss rates remain quite uncertain). 

Finally, we accounted for the effects of shock-generated X- 
ray emission in the model atmospheres. We adopted the standard 
Galactic luminosity ratios logLx/Lboi ~ -7 (Sana et al. 2006), 
with the bolometric luminosities determined in this study. 



3.2. Stellar and wind parameters 

To determine the stellar and wind parameters and surface abun- 
dances, we followed the same methodology as in previous pa- 
pers. We discuss here the main points only, and refer to Bouret 
et al. (2012) for a detailed description of our procedure and di- 
agnostics in more general cases. 

- Stellar luminosity: Since the COS and STIS spectra are 
flux-calibrated, we may use them together with optical and 
NIR photometry from Bonanos et al. (2010) (and refer- 
ences therein) to compute the luminosity and interstellar ex- 
tinction, by fitting the observed spectral energy distribution 
(SED) from the UV to the near-IR range. We used two cor- 
rections for reddening, one for the Galactic foreground us- 
ing the Cardelli et al. (1989) data and one for the SMC in- 
terstellar medium (Witt & Gordon 2000). The correspond- 
ing values are listed in Table 1. Both extinction curves are 
valid over the range 1 150 A to 30,000 A. The ratio of total- 
to-selective extinction was adopted to be Ry - 3.1 follow- 
ing Massey et al. (2009), while a distance modulus to the 
SMC of 18.91 ±0.02 (Harries et al. 2003) was used. Using 
Ry - 2.7 as in Evans et al. (2006) would induce changes 
of less than 5% in the flux, mostly in the UV, which is on 
the same order of magnitude or larger than the uncertainties 
caused by the photometric quantities or the distance modu- 
lus, and in any case would translate into uncertainties on the 
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luminosities of 0.04 dex or less. We found no indications 
for deviation from this standard Ry in our reddening correc- 
tion, which is sometimes suggested particularly in regions of 
strong star formation where Ry could be higher (e.g. in 30 
Dor, see Bestenlehner et al. 2011). We note, however, that 
Maiz Apellaniz & Rubio (2012) reported that four stars in 
a quiescent molecular cloud in the SMC presented signifi- 
cant variations in the extinction law, although all are located 
within a few pc of each other. The strongest variations were 
found in the 2175 A bump, and despite these UV variations, 
the SMC visible-NIR extinction law appears to be more uni- 
form, with no support found for high Ry. A somewhat more 
complex approach had to be adopted to derive the luminosity 
in three cases, which is discussed in Sect. 5.1. 

- Effective temperature and surface gravities: When opti- 
cal data were available, we derived an initial estimate of T e g 
with the He i /14471 to He n /14542 equivalent width ratio, 
as well as by fitting the line profiles. The T e s determination 
was then refined with other available He i and He n lines and 
listed in Table 2. Following Bouret et al. (2003), Hillier et al. 
(2003), and Heap et al. (2006), we also used the ionization 
ratio of the UV metal lines to estimate the effective temper- 
ature. This includes for instance the ratios from iron forest 
ions 2 , Crv/U169toCm/U176, or Nra AAl 183-1 185 and 
N hi XU748-1752 to N iv /H718 ratios. This is particularly 
important for the hottest stars of the sample, e.g., MPG 355, 
where the He i lines become very weak (see Rivero Gonzalez 
et al. 2012, for the use of nitrogen optical lines to constrain 
r ef f for the earliest O-stars). The C iv Al 169 to C in ^1 176 
ratio is especially useful because of its very good sensitivity 
to T e g and weak dependance on gravity, hence we used it as 
the primary T e ff diagnostic in the UV. However, in several 
cases, we found that effective temperature determinations 
based on the ionization ratio of helium lines are inconsistent 
with the ionization balance of carbon as predicted by CMF- 
GEN. In such cases, we found an offset toward higher r e jf by 
about 2000 K compared with those derived from the helium 
lines, and we were unable to achieve a good fit to these lines 
without degrading the agreement with the optical (helium) 
spectrum. The values listed in Table 2 refer to UV-based r e ff 
only for those stars for which no optical spectra were avail- 
able. For the other stars, the quoted T e g correspond to that 
providing the best fit over the whole FUV-optical spectral 
range. 

The Stark-broadened wings of Hy were the primary con- 
straint on log g, while the other Balmer lines were used as 
secondary indicators. The S/N of our spectra is sufficient to 
achieve a typical accuracy of about 0. 1 dex in log g. For stars 
with UV spectra only, we adopted a typical value for O-type 
dwarfs, log g- 4.0, because the UV range shows little sen- 
sitivity to changes in surface gravities (but see Heap et al. 
2006, for a discussion on using lines from different iron ions 
in the UV range). 

- Microturbulence: A constant microturbulent velocity of 15 
kms~' was adopted for calculating the atmospheric struc- 
ture. To calculate the emerging spectrum, we used a radi- 
ally dependent turbulent velocity (see also Hillier & Miller 
1998). The microturbulence in the photosphere (£, see Tab. 
2) is set such that it reproduces the strength of the Fe iv, Fe v, 
and S v lines in the COS spectra best. In the outer region of 
the wind, the turbulence was adopted to fit the shape and 



slope of the blue side of the absorption component of the 
C iv 1548-1551 P Cygni profile. 

- Surface abundance : We used SMC abundances from 
Hunter et al. (2007) for Mg and Si and from Venn (1999) for 
Fe. For all other elements but CNO, we adopted the standard 
solar abundances from Grevesse et al. (2007), scaled down 
by a factor of 5 (or equivalently by -0.7 dex). The strength 
of helium lines was used to constrain the helium abundance. 
We concentrated on deriving the carbon, nitrogen, and 
oxygen abundances. 

- Carbon abundance: In the UV range, we used C iv Al 169 
and C in Al 176 as prime indicators for carbon abundances, 
although these lines proved to be clearly sensitive to T s g and 
i>turb an d for the earliest stars of the sample showed some sen- 
sitivity to the wind density. For stars with later spectral types, 
C in A1247 shows a clear response to abundance variations. 
In the optical, we used C in i4070 and C iv ^5801-5806 
when possible. On the other hand, we did not include C in 
/15696 in our analysis. This line is sensitive to the adopted 
model atom and to UV C m lines at 386, 574, and 884 A 
(Martins & Hillier 2012). A fairly large number of energy 
levels is required for this line to form as an emission as ob- 
served in late-O stars, and Martins & Hillier (2012) called 
for extreme caution when interpreting carbon surface abun- 
dances derived from this line 3 . 

- Nitrogen abundance: In the UV, N in AAl 183-1 185 and 
N in /LU748-1752 are used as primary diagnostics. Together 
with N iv /U718, which is only marginally affected by the 
stellar winds even in the earliest stars, these lines are clearly 
photospheric and provide a very reliable estimate of N/H. 
N in /1/14634-42 and N iv /14058 are the strongest nitrogen 
lines in the optical but are often very weak in SMC stars, 
and were used as secondary diagnostics for nitrogen abun- 
dance. Their formation processes, including their sensitiv- 
ity to the background metallicity, nitrogen abundance, and 
wind-strength, have been discussed extensively by Rivero 
Gonzalez et al. (201 1, 2012). 

- Oxygen abundance: We used the O IV/LU338-1343 lines 
as a diagnostic for O/H in the UV, although they sometimes 
show a blue asymmetry that indicates a contribution from the 
wind for the earliest stars with stronger winds. We also used 
the O in AAl 150-1 154 A. These lines are on the shortest side 
of the UV spectra and the flux level there is more uncertain, 
however. O V/11371 was not used because it is weaker and 
is most sensitive to the mass-loss rate and clumping parame- 
ters. In the optical, we used O m /15592 when possible. 

- Wind parameters: The wind terminal velocities, v m , were 
estimated from the blueward extension of the absorption part 
of UV P Cygni profiles. The typical uncertainty in our de- 
termination of Voo is 100 kms -1 (depending on the maxi- 
mum microturbulent velocity we adopted). For stars with- 
our or with undeveloped P Cygni profiles, we scaled down 
the wind terminal velocities found in Galactic stars with the 
same spectral types. We started from values extracted from 
Table 1 in Kudritzki & Puis (2000), and applied a scaling v^ 
oc Z", with Z/Z = 1/5 and an exponent n - 0.13 (Leitherer 
et al. 1992). 

Mass-loss rates were derived from the analysis of UV 
P Cygni profiles. In contrast to the Galactic case (Martins 



2 Fe iv lines between 1550 and 1700 A, Fe v lines between 1300 and 
1550 A, and Fe vi lines below 1350 A. 



3 The influence of Fe metal lines found by Martins & Hillier (2012) 
should be less than for galactic O stars because the Fe metal abundance 
is down a factor of five 
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et al. 2012b), a unique value of M allows a good fit to both 
the UV lines and Ha, although there is a problem with N v 
,11238-1242 in some stars. The j3 exponent of the wind ve- 
locity law was derived from fitting the shape of the P Cygni 
profile. The clumping parameters, / and v c \ were derived in 
the UV domain following Bouret et al. (2003, 2005). Some 
photospheric lines in the optical presented some sensitivity to 
the adopted filling factor (and scaled M), For photospheric 
H i and He lines, for instance, this is essentially caused by a 
wind contribution that is weaker in clumped models, which 
produces deeper absorption than smooth-wind models. 
A different procedure had to be adopted for stars without 
wind signatures. Starting from mass-loss rates of Galactic 
stars of same spectral type, scaled to Zsmc = 1 /5 Z G with an 
exponent m - 0.83 (Mokiem et al. 2007) such that M cc Z m , 
these values were decreased until no P Cygni profiles were 
predicted for the aforementioned lines. As a consequence, 
and because we derived the mass-loss rates using homoge- 
neous models, the quoted M must be considered as upper 
limits. 

3.3. Rotation and macro-turbulence 

To compare the synthetic spectra to observations, we first con- 
volved the theoretical spectra with the appropriate instrumental 
profiles. Then, we convolved our synthetic spectra with rota- 
tional profiles to account for the projected rotational velocity 
(v sin i) of the star. In the UV range, the v sin i is constrained 
by matching relatively weak iron lines as well as C in A 1176, 
O iV/Ml 338- 1343, and S v ,11502. Helium lines were used in the 
optical. The C m Al 176 multiplet is especially useful because 
of its six close line components. Many spectra show partial or 
fully resolved multiplet components, indicating moderate (v sin i 
< 100 kms -1 ) to slow (usin i « 20 kms -1 ) rotational velocities, 
with a typical precision of 10 km s . A standard limb-darkening 
law was adopted for the convolutions with rotational profiles; a 
discussion about the impact of limb-darkening on some spectral 
lines of fast rotators is provided in Hillier et al. (2012). In most 
cases, rotational rates for the target stars had been measured in 
Penny & Gies (2009) from a cross-correlation function (CCF) 
analysis performed on FUSE spectra. The usin i we derived 
from a direct comparison to observed photospheric profiles are 
quoted in Table 2. They agree very well with those from Penny 
& Gies (2009), aand those by Mokiem et al. (2006) for the ten 
stars in common. 

Other broadening mechanisms are nowadays routinely in- 
cluded to account for (isotropic) macro-turbulence, i.e. symmet- 
rically distributed, deep-seated and highly supersonic velocity 
fields. Simon-Diaz et al. (2010) showed that a radial-tangential 
formulation of macro-turbulence fits the data better than a pure 
isotropic Gaussian. We used the Fourier transform method de- 
scribed in Simon-Diaz & Herrero (2007, and references therein), 
applied to wind-free lines in optical spectra to assess the rel- 
ative contribution of macro-turbulence and rotation to the ob- 
served broadening. Stars without optical spectra or with only 
low-resolution (i.e., 2dF) optical spectra were not included in 
this step. For all other stars, we found that usin i from line fit- 
ting and/or CCF analysis (Penny & Gies 2009) on one hand, and 
the values obtained from the Fourier transform analysis on the 
other hand, agree within the error bars. This agreement suggests 
a minor contribution from macro-turbulence. Since the respec- 
tive contribution of both broadening terms is quadratic, macro- 
turbulent velocity fields from a few kms -1 to several tens of 
kms may still be present, however. Aerts et al. (2009) empha- 



sized that using the Gaussian formulation of macro-turbulence 
to fit line profiles may lead to significant under estimation of the 
actual rotational velocity and warned that appropriate expression 
for the pulsational velocities need to be used instead. 

Therefore, the v sin i values quoted in Table 2 should be re- 
garded as lower limits, but we do not expect that this uncer- 
tainty significantly alters our conclusions for the surface abun- 
dances. First, the rotational term is always dominant with re- 
spect to macro-turbulence. Second, we note that the stars stud- 
ied here are dwarfs, and macroturbulent velocities seem to be 
significantly lower in dwarfs than in O supergiants (Markova 
et al. 2011, and references therein). Third, rotation and macro- 
turbulence in the adopted description are not expected to affect 
line equivalent widths, only the line profiles. 

4. Results 

The stellar and wind parameters we derived for the whole sample 
are gathered in Table 2. The effective gravity derived from the 
analysis were corrected for the effect of centrifugal forces caused 
by rotation (see also Sect. 5.4), following the approach outlined 
in Repolust et al. (2004). These log g c values were used to derive 
the spectroscopic masses. 

Chemical abundances by number relative to hydrogen are 
listed in Table 3 together with N/C and N/O abundance ratios 
and the adopted baseline for the abundance analysis. 

Given the complexity of the fitting process over the UV-to- 
optical spectral range, we elected not to derive formal statistical 
uncertainties, nor did we try to estimate the correlation of errors 
between parameters. Instead, we varied the different parameters 
until we obtained the solution that provided the best fit "by eye" 
to the observations, with more weight being given to primary 
diagnostics (see Sect. 3.2, and Bouret et al. 2012). Then, for 
a more quantitative estimate we varied each parameter indepen- 
dently around this solution and computed the residual (i.e., the 
difference) between the observed and the synthetic spectrum. We 
imposed the condition that this residual must remain within 10% 
and adopted the range of values that fulfilled this condition as 
error bars on each parameter. The stellar and wind parameters in 
Table 2 are obtained with this procedure. 

As for the uncertainty in the CNO abundance measurements, 
error bars are based on the fit quality of the lines listed above 
(Sect 3.2), as estimated from a x 1 procedure as presented for 
nitrogen in magnetic O-type stars by Martins et al. (2012a). Al- 
though this process is not ideal because it neglects correlations 
between parameters, it provides reasonable representative error 
bars. 

Appendix B provides a summary of our conclusions for each 
individual object, with best-model fits. Fig. 1 presents a selec- 
tion of UV spectra for stars representative of each spectral type 
in the sample, together with the corresponding best-fit models. 
Figure 1 shows that the complete observational sequence can 
be quantitatively reproduced by the models. Important photo- 
spheric features such as Fe iv to Fe vi ion lines or CNO lines 
are reproduced and strongly vary as functions of T e g, u tur b, and 
abundances. Wind profiles are clearly defined for stars earlier 
than 06 V. Note the steep transition from strong P Cygni pro- 
file at C iv /LU548-1550 to the characteristic weak-wind profile, 
occuring at spectral type 06.5 V and later. As seen from Fig. 
1, only a small fraction of stars in the present sample exhibit 
conspicuous wind signatures, which allows a good determina- 
tion of quantities such as the mass-loss rates, terminal velocities, 
and clumping factor (see also Table 2). We investigated these 
wind quantities (e.g. the wind-momentum-luminosity-relation) 
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Fig. 1. UV spectra for a selection of stars for each spectral type presented in this study, from earlier to later spectral types. Observations (black 
lines) are shown with best-fit model spectra (red lines). The observational sequence is very nicely reproduced by the models (Table 2). Note the 
strong transition in wind profiles (e.g. C iv /U1548, 1550) around spectral type 06 V. 



but too few stars provide reliable values for a fruitful discussion. 
For this reason, we chose not to discuss this in the present pa- 
per. We are currently investigating the giants and supergiants 
observed with COS/STIS during our program, and we postpone 
a complete discussion to a forthcoming, dedicated paper on wind 
properties at low metallicity. 

The qualitative behavior of the lines used for abundance de- 
termination in the UV is illustrated in Fig. 2. We selected stars 
with very similar fundamental parameters to show the abundance 
differences. For this comparison, the spectra were corrected 
for reddening (Table 1) and the presence of interstellar Lye ab- 
sorption was taken into account for the plots (mostly important 
for the spectral region where C iv A 1169, Cm All! 6 and N in 
AAl 183-1 1185 are found). The bottom panel of Fig. 2 shows 
AzV 243 and MPG 368, two stars with very similar parameters, 
including usin i and micro-turbulence (see Table 2), as shown 
by the similar S v A 1502 lines in the two stellar spectra. Car- 
bon and oxygen lines (C iv Al 169, C m A\ 176 and O iv /LU338- 
1 340) are very similar and indicate that carbon and oxygen abun- 
dances are the same in both stars. The nitrogen abundance on 
the other hand is slightly higher in MPG 368, as indicated by 
the stronger lines (marginally in the case of N m AAl 183-1 185 
but more clearly for N iv /11718, not shown here). The middle 
panel shows AzV 243 and AzV 446. The two stars have also 



very similar parameters as well, but for u sin i and u tu ,b, which 
are significantly higher in AzV 243. The CNO lines, although 
they look different in these stars, are best reproduced with simi- 
lar abundances. Finally, the top panel of Fig. 2 compares MPG 
113 and AzV 446, where the v sin i values are similar. The effect 
of the higher u tulb in MPG 113 is clearly visible in S v A15Q2. 
The CNO line profiles suggest that the surface abundance are 
different in these two stars: stronger carbon and lower nitrogen 
in MPG 113. Although O iV/Lll 338- 1343 are sensitive to v turb , 
the amplitude of the difference between the lines in the two stars 
suggests that oxygen abundance is higher in MPG 113 as well. 
We stress that the spectroscopic analysis did not indicate the 
need for oxygen abundances to differ from the baseline adopted 
for reference (Kurt & Dufour 1998) for about half of the sample 
stars (12 out of 23). 



5. Discussion 

5. 1 . Binaries and binary fraction 

In Sect. 3.2, we presented the method we used to constrain the 
stellar luminosity in the models. This method allowed us to 
reach very satisfactory fits of the global SEDs. Some remain- 
ing discrepancy between models and the flux-calibrated COS 
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Fig. 2. UV spectra for stars for similar fundamental parameters (Table 2). The profiles of C iv A 1169, C m /U176, and Nm AAl 183-1 1185 are 
shown in the left panels. Profiles of O iv /LU338-1340 and S v A1502 are presented in the right panels. 



spectra is nevertheless observed at wavelengths shorter than the 
Lya line. This discrepancy can be alleviated, to some extent, by 
adding an extra absorption caused by interstellar Lya. Fitting 
the wings of the observed Lya component on the COS spectra 
allowed us to derive an estimate of the neutral hydrogen column 
density on each line of sight. An additional problem is the un- 
certain extinction law in the SMC at these short wavelengths. 

For AzV 326, AzV 189 and AzV 148, the luminosities that 
were derived from the photometry, log L/L Q ^ 5.15, are 0.3 dex 
higher than expected for their spectral types 09V and 08.5V 
compared with stars of the same spectral type, either in the 
Galaxy (e.g. Marcolino et al. 2009) or in the SMC (e.g. Mok- 
iem et al. 2006). When used in conjunction with the effective 
temperatures and the surface gravities obtained from the spec- 
troscopic analysis, we found spectroscopic masses much larger 
than masses derived for stars of similar spectral type and lumi- 
nosity class, for instance in our Galaxy (Martins et al. 2005b; 
Marcolino et al. 2009). We stress that we were unable to find any 
mention of photometric and spectroscopic variability that would 
cause the observed shift in luminosities. A rather straightfor- 
ward explanation of these high luminosities is that these stars are 
actually binaries and that both the UV flux in the COS spectra 
and the optical/NIR photometry show contributions of the two 
components of the binary system. It is indeed not surprising to 
find some binaries in our sample, multiplicity being a common, 
almost ubiquitous feature among massive stars (Garcia & Mer- 
milliod 2001; Preibisch et al. 1999; Sana et al. 201 1). 

Furthermore, a detailed check of the Galactic O-type stars 
catalog by Maiz-Apellaniz et al. (2004) revealed that a signif- 
icant proportion of late-type binaries are found to be parts of 
couples of stars with the same spectral type (e.g., 08V + 08V, 
and 09V + 09V), hence each component contributes about half 



the observed flux. If we assume a similar behavior in the SMC, 
we could then propose the following correction. We divided the 
observed UV flux and fluxes from UB VIJHK photometry by two 
on one hand, and adopted a similar reduction of the luminosity in 
the model atmospheres. Since this simple procedure maintains 
an excellent fit of the overall SEDs (from the FUV to the NIR, 
see e.g. Fig. 3), it reliably supports the assumption that both 
stars have similar effective temperatures. This approach resulted 
in a reduction of the stellar radii, yielding spectroscopic masses 
that became then consistent with masses expected for late-type O 
stars. Moreover, the luminosities are also compatible with those 
expected for 08V - 09V Galactic stars (Martins et al. 2005a). 
The luminosities and masses listed in Table 2 are the values ob- 
tained with this approach. We emphasize that we do not argue 
that late-O stars in binary systems are systematically in equal 
mass systems, but simply that the three late O-type stars of our 
sample that show abnormal luminosities are more easily inter- 
preted in terms of such a scenario. Instead, the recent work by 
Sana et al. (2012a) for instance shows that low-mass stars might 
actually dominate the companion mass distribution. Reaching a 
more general conclusion on this specific question is beyond the 
scope of this paper and the sample it is based on. 

Three other stars, MPG 324, MPG 368, and MPG 487 have 
been claimed to be (potential) binaries. Evans et al. (2006) re- 
ported small radial velocity shifts in some lines in the spectra of 
MPG 324 (NGC 346-007) of about 20 - 30 km s 1 , suggesting 
it may be a single-lined binary. Massey et al. (2009) reported a 
velocity difference of 50 km s between the He i and He n lines 
of MPG 368, which they concluded is a strong indication of bi- 
narity. We note that a velocity difference of 70 km s _1 between 
the UV and optical spectra was also found by Heap et al. (2006) 
for this star. At the COS spectral resolution (=s 0.1 A), any corn- 
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Table 2. Fundamental stellar parameters. 
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-9.22 


0.8 




1950 


10 


50 


20.1+4.5 


NGC346-046 


07Vn 


39.0 


4.22 


4.81 


-9.22 


0.8 




1950 


10 


300 


19.0±4.3 


23.7 + *| 
23.2!| 


NGC346-031 


08Vz 


37.2 


4.00 


4.95 


-9.22 


0.8 




1540 


5 


25 


18.6+4.2 


AzV 461 
MPG 299 


08V 
08Vn 


37.1 
36.3 


4.05 
4.25 


5.00 
4.64 


-9.00 

-8.52 


0.8 
0.8 




1540 
1540 


10 

5 


200 
360 


23.8±5.4 
18.3+4.1 


23.8 + " 
19.5+lJ 

25 7 + " 

215^ 


3 6 +a6 

L5f 

4 3+8:? 

A 9+0.6 

-Q-5 
4 3 +(w 

-Q-6 
6.2 +lil * 

47+K6 

6.1 +0 ^ 
6.1 +0 ^ 
6 5 +0 ^ 

" -0.3 


MPG 487 a 


08V 


35.8 


4.00 


5.12 


-8.52 


0.8 




1540 


2 


20 


31.7+7.1 


AzV 267 c 


08V 


35.7 


4.04 


4.90 


-8.60 


0.8 


0.10 


1250 


10 


220 


21.9+5.0 


AzV 468 c 


08.5V 


34.7 


4.00 


4.76 


-9.15 


0.8 




1540 


5 


50 


15.7+3.5 


19.3 + |? 
19.0*|| 


AzV 148 a - c 


08.5V 


32.3 


4.00 


4.84 


-8.70 


0.8 




1540 


12 


60 


25.3+5.7 


MPG 682 


09V 


34.8 


4.10 


4.89 


-9.05 


0.8 




1250 


5 


40 


26.6±6.0 


20.8 + [5 
18.8 + ; : ° 
18.8*| : | 

21 + ^ 6 


AzV 326 a 


09V 


32.4 


4.03 


4.81 


-9.15 


0.8 




1250 


5 


200 


25.3±5.7 


AzV 189 a 
MPG 012 


09V 

BOIV 


32.3 
31.0 


4.02 
3.65 


4.81 

4.93 


-9.22 
-9.30 


0.8 
0.8 




1250 
1250 


5 
5 


150 
60 


24.8±5.6 
17.1+3.9 



Notes. (a) Stars whose photometry was interpreted as indicating binarity (cf. Sect 5.1) <b) Stars showing radial velocity variations (cf. Sect 5.1) 
(c) Stars whose log g - 4.0 was adopted, because we lacked of optical spectra. 

M are given in M yr~'. Evolutionnary masses (M cvo i) were derived from the tracks by Brott et al. (2011a) (see text for more comments). 
Uncertainties on T eff are ± 1500 K, ± 0.1 dex for logg (except AzV 446, AzV 429, AzV 267, AzV 468 and AzV 148, ± 0.2 dex), ± 0.1 dex on 
log L/L . As for the wind quantities, an uncertainty of ± 0.04 dex was estimated for stars with P Cygni profiles, while v^ was measured within ± 
100 km s~' . For stars without P Cygni profiles, see Sect. 3.2 for the procedure we adopted. 



1.0x10" 



2.2x10" 



4.6x10" 



1.0x10" 




0.55 
Wavelength (,u.m) 



3.0 



Fig. 3. Spectral energy distribution of AzV 189, from FUV (COS spectrum) to NIR (data from Bonanos et al. 2010). The fit to the SED with a 
model with a luminosity derived with the procedure outlined in the text, to take into account the potential binarity of this system, is shown in red. 
We used a distance modulus DM = 18.9 and E(B-V) Gal = 0.037 and E(B-V) SMC =0.1 
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ponent with radial velocity difference above 50 kms" 1 should 
be detectable, but we found no conspicuous evidence for SB2 
binaries in our FUV spectra. 

Bouret et al. (2003) argued that the weak UV lines of MPG 
487 might be diluted by one (or several) additional stars in the 
slit. From the present modeling, this conclusion was somewhat 
re-inforced by the failure of fitting the observational SED, from 
the FUV to the NIR, without down-scaling the optical-NIR pho- 
tometric fluxes, by 10% to 20% over the parameter space we ex- 
plored. The best fit we obtained (see Table 2 for the parameters) 
implied a reduction of 10% of these fluxes, but we did not need 
to change the UV flux accordingly. This may be the result of a 
companion with a significantly later spectral type and a contribu- 
tion to the total luminosity of only about 10% of that measured 
for MPG 487. Because it is cooler, it would mostly contribute 
to the optical and NIR bands and the lower luminosity would 
make the change in the photometry within the observed range. 
We argue therefore that this suggests that MPG 487 may be in a 
binary system. Of course, effects on the luminosity or the SED 
may result from the presence of several objects in the instrumen- 
tal aperture, and cannot be unambiguously tied to physical bina- 
ries as for systems identified through radial velocity variations. 
Furthermore, the fact that we are able to match luminosities and 
SEDs of the other stars with single-star models, does not exclude 
the possibility that these stars have cooler and fainter compan- 
ions that remain undetected in the SED. Our study is obviously 
not based on a complete and homogeneous population of O-type 
stars in the SMC, because it mixes field stars and stars in the 
open cluster NGC 346. It is biased towards the detection of rel- 
atively massive companions, that would contribute most to the 
total observed flux (AzV 326, AzV 189, AzV 148, MPG 487) 
or induce radial velocity variations of several tens of kms" 1 on 
early-type objects (MPG 324 and MPG 368). To summarize, 
although we did not find additional hints for binaries from the 
modeling analysis, it is not unlikely that there are more multiple 
objects in our sample than the six stars discussed above. Our 
estimate of the binary fraction should therefore be viewed as a 
lower limit of the actual fraction of multiple systems in the SMC 
that includes O-type stars. 

With 6 binaries in a sample of 23 stars, the binary frac- 
tion would be 26%. This is consistent with recent results on 
the fraction of binaries among massive stars in the Magellanic 
Clouds. The first VLT-FLAMES survey of massive stars re- 
ported in Evans et al. (2006) found a lower limit to the binary 
fraction for O stars in NGC 346 of 21% (the total fraction for 
O-type and early-B-type stars being 26%), while the investiga- 
tion of multiplicity in LMC massive stars by Sana et al. (2012a) 
found that the observed binary fraction is 35% but the intrinsic 
fraction is as high as 5 1 %, after correcting for observational ef- 
fects. 

5.2. H-R diagram, evolutionary status, and ages 

Using the effective temperature and luminosity from the spectro- 
scopic analysis, we built the H-R diagram shown in Fig. 4. For 
this diagram, we also used evolutionary tracks and isochrones 
from Brott et al. (2011a) for initial rotational velocities « 180 
kms" 1 . The following considerations lead us to select this later 
value: First, the average wsin i of our sample is « 115 kms"' 
(for a standard deviation of 1 kms" 1 ). We recall that as dis- 
cussed in Sect. 2, part of our sample (stars with STIS spectra) 
was selected from slow rotators and therefore the average v sin i 
should be regarded as biased towards lower values. Second, from 
their subsample of (21) unevolved OB stars in the SMC, Mok- 



iem et al. (2006) concluded that the average v rot of their sample 
was 150-180 kms" 1 . These two points made us choose use the 
u rot of the Brott et al. (2011a) models which is most consistent 
with the mean velocity from Mokiem et al. (2006), that is, u rot 
= 180 kms" 1 , although it must be kept in mind that this u rol 
for evolutionary models is the initial rotational velocity, without 
projection effects. The error bars on the position of stars in the 
H-R diagram reflect the uncertainties on the effective tempera- 
ture and luminosities from the analysis above. Fig. 4 shows that 
the stars have ages ranging from younger than 1 Myr to 7 Myr. 
Their evolutionary masses range from 19 M Q to 42 M Q with the 
exception of MPG 355 which is apparently more massive than 
60 M e , consistent with its early spectral type and correspond- 
ing high effective temperature and luminosity. Our sample con- 
sists of field stars as well as stars belonging to the giant H II 
region and to the open cluster NGC 346. Overall, the age span 
of the two groups overlap, but we note that no stars younger than 
1 .5 Myr are found in the field sample, suggesting a lack of recent 
massive star formation in the SMC field. 

Of the twelve stars in NGC 346, two are on or very close to 
the ZAMS. A third, rapidly rotating (see below) object is found 
between the isochrones one and two Myr (Table 2). A second 
group of five stars are between two and four Myr old, while two 
objects are younger than 5 Myr. The age spread is very consis- 
tent with the previous study of Mokiem et al. (2006). Finally, 
the coolest star of our sample, namely MPG 012, is older than 
6 Myr. Our analysis, which includes the FLAMES and UVES 
spectra of MPG 012, indicates that this star is the oldest of the 
sample, significantly older than all analyzed stars in NGC 346 
(see also Bouret et al. 2003). The evolved status of MPG 012 is 
also supported by its lower surface gravity (as also suggested by 
the luminosity class IV ascribed to this star) and a high enhance- 
ment in surface nitrogen. These findings support the argument 
of Walborn & Howarth (2000) that MPG 012 is most likely not 
a coeval member of the cluster, which the based on its spatial lo- 
cation outside the nebula, a somewhat discrepant radial velocity, 
and the derived stellar parameters. 

MPG 355 is the hottest and most luminous star of our sam- 
ple. The optical spectra are devoid of He i lines and we had to 
rely on alternative diagnostics to constrain the effective temper- 
ature (but see Massey et al. 2009, for a mention to a different 
spectrum, that potentially shows a very weak He i /14471 line). 
Our parameters for this star agree very well with those from the 
recent analysis of O-stars of the earliest spectral types by Rivero 
Gonzalez et al. (2012). Bouret et al. (2003) argued that the posi- 
tion of MPG 355 close to the ZAMS could indicate chemically 
homogeneous evolution due to rotation close to break-up veloc- 
ity (Maeder 1987). Testing different models with high initial 
masses (50 and 60 M G ) and high initial rotation rates (above 400 
kms" 1 ), we found that the isochrone 3 Myr for models with ini- 
tial u ro , = 400 km s" 1 comes closest to the location of MPG 355 
in the H-R diagram (see Fig. 5). Still, MPG 355 is above this 
isochrone and above the 60 M G , u rot = 400 kms" 1 evolutionary 
tracks, and the main conclusion we can reach at this point is that 
MPG 355 is a very massive object, probably more massive than 
70 M Q , but its actual evolutionary status is more uncertain. Ded- 
icated evolutionary models are needed to reach a better quanti- 
tative conclusion regarding its status. We note that a chemically 
homogeneous evolution was concluded to be unlikely in Rivero 
Gonzalez et al. (2012) because of the low usin i measured for 
this star. The analysis of the surface abundances support a non- 
standard evolutionary status for this object, however (see Sect. 
5.3). 
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Fig. 4. H-R diagram of the sample stars. Evolutionary tracks and isochrones are taken from Brott et al. (201 la) for an initial rotation rate of 180 
kms -1 . Evolutionary tracks are annotated in solar masses, and isochrones in Myr from 1 to7 Myr. Green circles indicate field stars, while squares 
indicate stars in NGC 346. Blue squares show normal luminosity class V stars, while black squares show Vz-class objects (cf. Sect. 5.2.2). Red 
squares show the position of the two stars that rotate faster than 300 kms -1 . 



Overall, there is a trend for the most massive stars in the sam- 
ple to be younger than less massive stars. However, it is not clear 
how much of this finding is true evidence for a real age sequence 
with mass or instead is due to a selection effect of the sample, 
because we selected stars to be primarily dwarfs. Indeed, more 
massive stars have shorter lifetimes and are more quickly identi- 
fied as evolved, and hence not as dwarfs. In other words, above 
a certain mass, any massive star identified as a dwarf must be 
somewhat younger than less massive counterparts. For instance, 
an O star more massive than 40 M has already moved a long 
way from the ZAMS at an age of 3 Myr (see Fig. 4) and any 
star more massive would appear as a giant or even supergiant, 
hence it would not have been selected in the sample. It could 
also be a bias from the spectroscopic analysis, as observational 
results found the same trend for Galactic O-type dwarfs in clus- 
ters (Martins et al. 2012b) or in Galactic supergiants (Bouret et 
al. 2012). 

If true, however, in the specific case of NGC 346 this trend 
would suggest a scenario in which the most massive stars of the 
cluster formed last, and after the had formed, they would quench 
subsequent star formation (as suggested originally in Heap et al. 
2006). This trend can also be seen in the spectroscopic study of 
NGC 346 by Mokiem et al. (2006) (see their Fig. 12), although 
not as clearly. When comparing the stars in common with our 
analysis and the work by Mokiem et al. (2006), we note that they 
obtained effective temperatures systematically higher than ours 



(1700 K on average). This systematic difference has two impor- 
tant consequences: i) their sample is globally shifted toward the 
ZAMS compared to ours, it) most stars lie on a relatively narrow 
band, especially in the lower part of the H-R diagram, which 
is much less sensitive to age because of the smaller separation 
between isochrones. 

An explanation for the observed trend might be that the 
tracks adopted to derive an age for higher and lower masses stars 
should not have the same rotation rate. If for some reason the 
initial rotation rate v ro , of more massive objects turned out to 
be higher than for lower-mass stars, tracks corresponding to this 
higher v rot should be used preferentially. It is now a strong con- 
clusion from stellar evolution models with rotation that the main- 
sequence lifetime increases as rotation is faster (e.g. Maeder & 
Meynet 2000; Brott et al. 2011a). Age determinations based on 
non-rotating models can be biased toward younger age, which in 
contrast means that for models with higher rotation rates, older 
ages would be derived. On the other hand, actual differences be- 
tween age estimates from non-rotating and rotating models with 
moderate rotation should remain smaller than the error bars re- 
sulting from the uncertainties on temperature and luminosities 
for the relatively unevolved objects discussed here. In the spe- 
cific case of NGC 346, Mokiem et al. (2006) investigated the im- 
pact of using non-rotating versus rotating Geneva models (with 
initial v rot - 300 kms~' from Meynet & Maeder (2005)) to de- 
rive ages of the cluster stars. For stars on the main-sequence, 
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the authors estimated that the age difference may increase up to 
0.5 Myr (see their Fig. 13 and 14). We concur with this con- 
clusion using models without rotation and with initial rotational 
velocities of » 200 km s from Brott et al. (201 la). 



5.2.1. Fast rotation 

The fastest-rotating two stars of the sample, NGC 346-046 and 
MPG 299, lie close to or on the ZAMS in the lower part of the 
H-R diagram in Fig. 4. The derived evolutionary masses are 

23.v\i M ° and 19 - 5 -!o M o and a § es of °- 3 M y r and L5 M y r ' 

respectively, for models with v rol w 180kms _1 . It is well-known 
since the work of von Zeipel (1924) (see also Maeder 1999) that 
a fast-rotating star has stronger radiative fluxes at the pole than at 
the equator. Therefore, the observed luminosity, effective tem- 
perature, and also the effective gravity of a star depend on its 
initial angular velocity as well as on the inclination angle. In 
turn, its location on the H-R diagram depends on the same quan- 
tities, hence its age and mass as determined by fitting evolution- 
ary tracks or isochrones. This holds especially for stars that are 
located near the zero-age main sequence. A young, unevolved, 
slowly rotating star may have the same effective temperature and 
luminosity as a less massive, fast-rotating, evolved star. 

Brott et al. (2011a) showed that for fast-rotating stars (with 
initial v ro , = 400 kms 1 or higher), the location on the H-R dia- 
gram stays closer to the ZAMS than for slower rotators of same 
age. For this extreme rotation, mixing is so efficient that the 
whole internal structure and composition of a star is modified 
and may lead to homogeneous evolution in which the star be- 
comes brighter and hotter and evolves up- and blueward directly 
in the H-R diagram. The onset of a blueward bifurcation of the 
evolutionary tracks caused by chemically homogeneous evolu- 
tion arises around 18 M Q at SMC metallicity for stars with v rot 
» 550 kms~' or higher. This strongly suggests that evolutionary 
tracks with such high rotations rates are better suited to estimat- 
ing the evolutionary status of NGC 346-046 and MPG 299. In 
Fig. 5, we show that the location of NGC 346-046 and MPG 299 
on the H-R diagram can be accounted for by various isochrones, 
with combinations of ages and rotation rates such as (1 Myr, v rot 
« 400 kms -1 ) or (5 Myr, v roi ~ 550 kms~'). In other words, 
using isochrones with lower initial v ro , leads to a significant un- 
derestimation of the age of NGC 346-046 and MPG 299. 

It is also argued that rotation induces an enhancement of 
mass loss, which should affect the ultimate evolution of a star 
(Owocki et al. 1996; Maeder 1999; Maeder & Meynet 2000). 
However, we calculated that a contrast of 10% at most should 
be observed for the mass-loss rates at the equator and at the pole 
(see Maeder & Meynet 2000, for a discussion), for the radii, 
masses, and osin i we derive for NGC 346-046 and MPG 299, 
indicating that the actual increase of mass loss by rotation is very 
marginal for these two stars. The net effect of this an increase on 
the observed spectra is also expected to be negligible, which is 
also supported by the very good quality of the fits we obtained 
for the UV and optical spectra of both stars, indicating no sig- 
nificant deviation from spherical symmetry of their winds. The 
mass-loss rates we derive for NGC 346-046 and MPG 299 are 
similar to those of stars with same spectral type and much lower 
usin i, which belong to the range of spectral types where the 
weak-wind phenomenon is at play. In any case, a change by a 
factor of a few is probably not enough to be detected even in the 
UV spectra of these stars. This result implies that these stars lost 
only a marginal fraction of their initial angular momentum dur- 
ing the course of their evolution. Therefore, their rotation rates 
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Fig. 5. H-R diagram showing evolutionary tracks (red dashed line) 
with 60 M e and 20 M Q for an initial rotation rate of as 400 km s _1 (Brott 
et al. 2011a). The ZAMS is indicated by the black full line and the 
isochrone 5 Myr is indicated by the black, dot-dashed line. The green 
and blue line show the isochrones 1 Myr and 3 Myr, respectively. Blue 
filled square indicates the position of MPG 355, while the red squares 
show the positions of the two stars that rotate faster than 300 km s" 1 . 
The black square indicates the position of MPG 523. 



should not have decreased much and the present measured v sin i 
should reflect the initial values well. 

Mokiem et al. (2006) found that four SMC O dwarf stars 
close to or on the ZAMS had an helium-rich surface, two of 
them being the fast rotators discussed in this Section. Since 
helium enrichment is expected for fast rotators, even early in 
the stellar evolution (Meynet & Maeder 2000), Mokiem et al. 
(2006) concluded that the two stars with high helium surface 
abundance (Yne - 0.16) and low usin i were actually fast rota- 
tors seen pole-on. The four stars are included in our own work 
(MPG 113, NGC 346-031, NGC 346-046, and MPG 299). Us- 
ing the surface helium abundance as a measure of stellar age 
for a given mass for of stars with very high initial v ro u Mokiem 
et al. furthermore derived lower limits on the ages of these stars 
that were significantly higher than for non-rotating (or slowly- 
rotating) stars. The authors concluded that these stars probably 
follow a chemically homogeneous evolution. 

We used the same FLAMES spectra and supplemented them 
with UVES spectra for the first three stars. However, we found 
no evidence of a helium abundance above the initial value, ex- 
cept for MPG 299. In this case only, we agree with Mokiem 
et al. (2006) that an enrichment up to Y He = 0.13 provides a bet- 
ter match to the optical spectrum. Given the high v sin i of this 
star, this is indeed suggestive of chemically homogeneous evolu- 
tion. The conclusion of Mokiem and collaborators that MPG 113 
and NGC 346-03 1 are fast rotators cannot be supported from our 
analysis. 



5.2.2. The luminosity class Vz 

Three stars in our sample, MPG 113, MPG 523, and NGC 346- 
031, are classified as Vz stars (Evans et al. 2006; Walborn & 
Howarth 2000). This luminosity class was introduced by Parker 
et al. (1992) and Walborn & Blades (1997) for stars whose prop- 
erties are expected to be characteristic of extreme youth. The 
defining spectroscopic criterion for this class is that at every 
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spectral type, He n /14686 absorption is stronger than any other 
He line. More specifically, this would be the case of He n /14541 
at early-O types, but at type 07, the latter is equal to He i /14471, 
so that He n /14686 is stronger than both in a Vz spectrum. At 
later-O types, however, He n /14541 weakens more rapidly with 
advancing type than He n /14686 in normal class-V spectra, such 
that, a late-O Vz spectrum must have He n /14686 stronger than 
He i /14471. As stated in an illuminating discussion by Walborn 
(2009), this a behavior very likely results from the fact that for 
objects in the Vz class, a deficit of emission filling in He H/14686 
is observed, in contrast to normal class-V stars. This deficit re- 
sults from an "inverse" Of effect (see e.g. Walborn 2001, for a 
thorough definition), thus implying that O Vz stars should be 
subluminous and less evolved. They are therefore expected to be 
close to or on the ZAMS. 

Of these three stars, only MPG 523 is found to be located 
on the ZAMS within the error bars; the other two stars are close 
to the 3 Myr isochrone (see also Table 2). In the previous study 
of Mokiem et al. (2006), MPG 523 was clearly to the left of 
the ZAMS because of a derived T e s higher by 3000 K and only 
slightly higher log L/L G (by 0.07 dex). The strong nebular emis- 
sion present in the FLAMES spectrum of this star can be inter- 
preted as another hint at the young age of this star. On the other 
hand, Fig. 5 shows that the location of MPG 523 in the H-R 
diagram is actually consistent with that of a 20 M G star, rotating 
near 400 km s that is either 1 Myr or 5 Myr old. This interpre- 
tation is also confirmed by the spectroscopic analysis of MPG 
523, which indicates surface abundances in stark contradiction 
with the expectation for a young, slowly rotating star. From its 
abundance patterns, MPG 523 would be more compatible with 
being a more evolved, fast-rotating object. 

The properties of MPG 113 are also very non-standard. First, 
we did not find this star to be close to or on the ZAMS, in con- 
tradiction with Mokiem et al. (2006). The observed shift in the 
H-R diagram is again caused by an effective temperature about 
3000 K higher than our value, while luminosities agree to within 
0.05 dex. We recall that the helium surface abundance of MPG 
113 corresponds to that of a normal main-sequence 06V star, 
with no indication for enrichment {Yn e - 0.16 is firmly ruled 
out). Second, the carbon and nitrogen contents of MPG 1 13 are 
higher than the values measured in H II regions of the SMC (Du- 
four & Kurt 1998), resulting in an N/C ratio that is significantly 
higher than the standard ratio measured in H II regions (in agree- 
ment with Bouret et al. 2003). This does not support the claim 
by Walborn & Howarth (2000) that the N/C abundance ratio in 
MPG 113 is atypically low to explain its spectral classification 
as an OC6 star. We also note that another "normal" 06.5 V star 
has even higher carbon abundance but is also richer in nitrogen, 
yielding the same N/C ratio, although it is classified as an OC6 
star. It seems more likely that this classification is more complex 
than a simple abundance criterion and involves several physi- 
cal parameters such as the effective temperature and the surface 
gravity. 

For the third star, NGC 346-031, although we found it again 
to be cooler than did Mokiem et al. (2006) for very close lumi- 
nosities (within 0.04 dex), we found that this star does not lie 
especially close to the ZAMS and that it is older than 2.5 Myr. 
Mokiem et al. (2006) argued that the Vz classification, from a 
strong He n /14686, might come from the high He abundance 
they derived (Yg e - 0.16). We needed no such overabundance in 
our own models to obtain very good fit for this star (see Fig. 6). 
Indeed, using Yn e as high as in Mokiem et al. (2006) degrades 
the quality of the fit significantly. On the other hand, we de- 
rived a mass-loss rate that is two orders of magnitude lower than 
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Fig. 6. Wind lines of NGC 346-03 1 . A comparison of our best-fit 
model (blue line) with a model computed using the parameters provided 
in Mokiem et al. (2006) (red line). The mass-loss rates differ by almost 
two orders of magnitude. The blue side of Nv/il240 is affected by 
interstellar Ly a absorption. We also show the Ho- line, which proves to 
be insensitive to mass-loss over a wide range of values. 



that of Mokiem and collaborators, from the modeling of the UV 
spectrum (Fig. 6). The contribution of such a small M to a po- 
tential wind-emission that would come and fill in the He n /14686 
line is drastically reduced. We recall that for stars without con- 
spicuous PCygni profiles, we used homogeneous wind models, 
and that the M we quoted are upper limits above which models 
start producing wind signatures not seen in the data. We note 
that the three Vz stars have much smaller M in our UV based 
study than do those from Mokiem et al. (2006), which results in 
all cases from the extreme weakness (or absence) of the observed 
UV stellar-wind profiles. In the temperature range of interest for 
these stars, this might be the reason why the He n /14686 is ob- 
served to be stronger than in normal dwarf stars that experience 
stronger wind emission. For later-type stars, a lower temperature 
would act such that He i lines remain stronger than He n /14686 
even in cases where weak-wind emission is involved. 

The defining spectroscopic criterion for the Vz class seems 
to fail in establishing a one-to-one correspondence between an 
object and a young evolutionary status (but see Walborn et al. 
2013, for recent findings from VFTS survey) as inferred from 
the location in the H-R diagram. The bias toward lower T e g (than 
optical-based r e ff) to reproduce the observed ratio of UV carbon 
lines is an obvious cause for the shift of this location away from 
the ZAMS. The same trend was already noticed, with almost the 
same amplitude, for UV-based T e s derived with Tlusty in Heap 
et al. (2006). It is therefore unlikely that this is related to the 
physical assumptions we made (e.g., including a stellar wind, the 
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adopted micro-turbulence, the atomic data or other issues related 
to the sphericity of the atmosphere as considered in CMFGEN). 
The same trend was also found by Garcia & Bianchi (2004) from 
a UV-based analysis, although in a different context. In any case, 
our modeling suggests that more criteria than the strength of He 
lines should be taken into account to infer the true evolutionary 
status of a star, including the surface abundance patterns (see 
also Sect. 4. Brott et al. 201 la, for a discussion). 

5.3. Chemical evolution 

While rotation was claimed to be able to account for the early 
chemical enrichment of the stellar surface (e.g. Maeder & 
Meynet 2000), several observational studies have undermined 
this conclusion. Investigating the link between surface abun- 
dance patterns and rotation, either in the Galactic context or in 
the Magellanic Clouds, Hunter et al. (2007, 2008, 2009) con- 
cluded that up to 40% of their sample presented nitrogen surface 
abundances that could not be accounted for by rotational mixing 
in stellar evolution models. The stars challenging the link be- 
tween surface abundances and rotation fall into two categories: 
first, stars with strong enrichment but slow rotation, and second, 
fast rotators with a low (if any) surface enrichment. These stud- 
ies focused on nitrogen surface abundances of B-type stars up 
to the supergiants stage. In the core-hydrogen burning phase, 
slowly rotating nitrogen-rich stars are found for the lower metal- 
licity regime of the LMC and SMC (Hunter et al. 2009), while a 
relative deficit of such stars is observed in the Galactic case. 

This strongly suggests that the global chemical evolution of 
massive stars, as probed from their surface abundances, is more 
complex than assumed and is an intricate combination of sev- 
eral physical parameters, including mass, metallicity, age (evo- 
lutionary status) and rotation (Maeder et al. 2009). Martins 
et al. (2012b) showed that Galactic O stars in NGC 2244 and 
Mon OBI populated the entire nitrogen abundance - u sin i dia- 
gram, but nicely followed the predictions of evolutionary models 
in the nitrogen abundance - luminosity diagram, confirming the 
arguments of Maeder et al. (2009). Another obvious complica- 
tion is binarity. Using population synthesis models, Brott et al. 
(201 lb) concluded that the observed distribution of Hunter et al. 
(2008) could not be reproduced by single-star population synthe- 
sis. However, they also concluded that the nitrogen-rich slowly 
rotating stars are unlikely to result from close binary system evo- 
lution, with mass transfer and tidal spin-down. A fraction of this 
population would end up overpopulating the group of nitrogen- 
rich fast rotators, which is not observed. 



5.3.1. Helium abundances 

For massive stars in the hydrogen-burning phase of their evo- 
lution, helium and nitrogen are produced at the expense of car- 
bon first, as the star progresses to CN-equilibrium. The rela- 
tive increase of abundances is expected to be weaker for helium 
than for nitrogen, however, because the former is already abun- 
dant in the atmosphere while the latter is not. This is confirmed 
by stellar evolution models with rotation (Brott et al. 2011a), 
which predict that helium should remain fairly constant during 
the main-sequence phase and that helium transport to the sur- 
face is only substantial for the more massive and fast-rotating 
stars where mixing processes are more efficient. Quantitatively, 
the alteration of the helium surface abundance remains below 
0.2 dex during the H-burning phase for stars with M+ < 25 M e . 
Only fast-rotating models (u rot > 330 kms -1 ) for stars 30 M Q 



and above present modification of the helium surface content on 
the main-sequence larger than 0.25 dex. Our spectroscopic anal- 
ysis only found the need for helium enrichment in two stars (see 
Table 3), and for these two stars the enhancement was higher 
than expected. Of these two stars, one is indeed a confirmed fast 
rotator (MPG 299), the other a suspected binary (MPG 487). 
Overall, the qualitative and quantitative agreement of our mea- 
surements with theoretical predictions is therefore excellent. 

Since the enhancement of the surface helium content is low, 
a comparison of the predicted versus observed values is far from 
conclusive for such an early evolutionary phase. In principle, a 
better diagnostic is the N/He ratio, which quickly evolves from 
the initial value to a value corresponding to the CNO equilib- 
rium. Stars for which the equilibrium has not been reached will 
present intermediate values of N/He. Here again, however, the 
surface enrichment in helium is such that the N/He ration mer- 
erly witnesses the global surface enrichment in nitrogen, which 
we discuss in more detail in the next sections. 

5.3.2. Nitrogen abundances and comparisons with B stars 

Our results on the surface abundances of O stars in the SMC can 
be viewed as an extension of the work of Hunter et al. (2009) to- 
ward higher mass stars. In Fig. 7, we show the e(N) - v sin i dia- 
gram for our sample stars and the B stars (core hydrogen burning 
sub-sample) of Hunter et al. (2009). At first sight, the distribu- 
tion of O and B stars in this diagram is the same. This is espe- 
cially true for the low-velocity part of the diagram (u sin i < 100 
km s -1 ) where actual abundances have been measured in most B 
stars, not just upper limits (which is the case for fast-rotating B 
stars). For all masses, the highest value of e(7Y) seems to be about 
7.9. This suggests that the mechanisms responsible for the chem- 
ical enrichment of slowly rotating massive stars depends weakly 
on the mass. In other words, for any chemically enriched, slowly 
rotating massive star, the degree of enrichment is controlled by a 
parameter that is neither the mass nor the rotation rate. However 
we cannot exclude that the main-sequence B stars from Hunter 
et al. (2009) are older than the O stars of our sample. Since rota- 
tional mixing is expected to have a stronger effect toward core- 
hydrogen exhaustion, the observed match of B-type and O-type 
distributions could be spurious because it may be biased by the 
fraction of the main-sequence lifetime actually covered by each 
star (see e.g. Fig. 12 in Brott et al. 201 la). 

In the right part of Fig. 7 (above 100 km s~'), O and B stars 
also populate the same region. However, it is important to note 
that the values of e(N) for the B stars are mostly upper limits. 
The only safe conclusion one can draw is that O stars are at least 
as enriched as B stars. But one cannot exclude that O stars are 
substantially more nitrogen-rich. This result would be consistent 
with predictions of evolutionary tracks (Maeder & Meynet 2000; 
Brott et al. 201 la; Ekstrtim et al. 2012) and the results of Martins 
et al. (2012b), who showed that Galactic O-stars seem to have 
higher surface nitrogen abundance than B stars. 

In Fig. 8 we compare of the observed surface nitrogen abun- 
dance of O stars with predictions of evolutionary tracks. We 
chose the tracks of Brott et al. (2011a) since they are currently 
the only ones computed for different rotational velocities at the 
SMC metallicity. We divided our sample into two groups: stars 
below and above 27 M Q , the average mass of our sample. In the 
lower panel of Fig. 8, we show the low-mass part of our sam- 
ple. The group of slowly rotating N-rich stars is not reproduced 
by the tracks, since those tracks reach the low-velocity - high N 
regime only well beyond the end of the main sequence, repre- 
sented by the sudden left bifurcation in the tracks. This confirms 
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Table 3. He, C, N, and O abundances by number. 



Star 


Y He 


KC) 


e(N) 


6(0) 


N/C 


N/O 


MPG 355 


0.09 


7.39±0.15 


8.11+0.16 


8.00+0.18 


5.25 


1.29 


AzV177 


0.09 


7.17±0.16 


8.01+0.22 


7.90+0.30 


6.92 


1.29 


AzV 388 


0.09 


7.39+0.15 


7.96+0.20 


8.05+0.21 


3.71 


0.81 


MPG 324 


0.09 


7.00+0.15 


7.38+0.20 


8.00+0.24 


2.40 


0.24 


MPG 368 


0.09 


7.20+0.14 


7.70+0.20 


8.05+0.18 


3.16 


0.45 


MPG 113 


0.09 


7.52+0.16 


7.40+0.34 


8.30+0.12 


0.76 


0.13 


AzV 243 


0.09 


7.20+0.20 


7.50+0.34 


7.90+0.12 


2.00 


0.40 


AzV446 


0.09 


7.20+0.22 


7.48+0.32 


7.98+0.14 


2.04 


0.32 


MPG 356 


0.09 


7.69+0.17 


7.38+0.30 


7.98+0.08 


0.49 


0.25 


AzV 429 


0.09 


7.69+0.17 


7.08+0.29 


7.98+0.22 


0.25 


0.13 


MPG 523 


0.09 


6.38+0.15 


7.08+0.35 


7.40+0.10 


5.00 


0.48 


NGC346-046 


0.09 


6.68+0.14 


7.56+0.45 


7.98+0.37 


7.59 


0.38 


NGC346-031 


0.09 


7.69+0.10 


7.08+0.35 


7.98+0.08 


0.25 


0.13 


AzV 461 


0.09 


7.38+0.15 


7.38+0.26 


7.98+0.35 


1.00 


0.25 


MPG 299 


0.13 


7.69+0.15 


7.56+0.44 


8.10+0.33 


0.74 


0.29 


MPG 487 


0.15 


7.39+0.20 


6.60+0.26 


7.98+0.08 


0.16 


0.04 


AzV 267 


0.09 


7.30+0.14 


7.48+0.30 


7.98+0.30 


1.51 


0.31 


AzV 468 


0.09 


7.30+0.17 


6.30+0.25 


7.60+0.14 


0.10 


0.05 


AzV 148 


0.09 


7.69+0.15 


7.26+0.30 


7.98+0.10 


0.37 


0.19 


MPG 682 


0.09 


7.39+0.14 


6.78+0.25 


7.98+0.14 


0.25 


0.06 


AzV 326 


0.09 


7.08+0.14 


7.38+0.30 


7.98+0.24 


2.00 


0.25 


AzV 189 


0.09 


7.69+0.15 


7.38+0.30 


7.98+0.21 


0.49 


0.25 


MPG012 


0.09 


7.34+0.16 


7.70+0.21 


8.05+0.16 


2.29 


0.45 


SMC(HII) 


0.09 


7.37 


6.50 


7.98 


0.13 


0.03 



Notes. For comparison, we list the chemical composition determined by Kurt & Dufour (1998) for HII regions in the SMC as used in the 
evolutionary models by Brott et al. (201 la). By convention log e(X) = 12 + log [X/H]. 
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Fig. 7. Nitrogen abundance (in units of 12+log(N/H)) versus projected 
rotational velocity. The black squares indicate the B stars of Hunter 
et al. (2009). The red circles show the O stars analyzed in this study. 



the well-known results that rotational mixing does not explain 
the high N content of some of the slowly rotating massive stars. 
O stars rotating faster than 100 kms~' can be accounted for by 



evolutionary tracks provided the stars are all very close to the 
end of the main sequence. This situation is rather unlikely, and 
we must therefore conclude that a slightly stronger mixing in the 
models would be necessary to explain the range of observed N 
abundances. 

This conclusion applies even more strongly to the high-mass 
subsample displayed in the upper panel of Fig. 8. In the high- 
velocity part, there are only three stars, but they all fall well 
above the main-sequence part of the tracks. Either the evolution- 
ary models do not predict enough mixing, or the stars are already 
in a post-main-sequence phase. The latter possibility is rather 
unlikely given their observed properties. It is now well estab- 
lished that very massive stars can still be core-H burning objects 
and appear as Wolf-Rayet objects (Martins et al. 2008; Crowther 
et al. 2010). Here, one would need the exact opposite to ex- 
plain the strong N enrichment of massive SMC O stars: stars 
already in a post-core-H-burning phase, but with the appearance 
of main sequence stars. This plot thus favors the possibility that 
rotational mixing is not sufficient at high masses in the models 
of Brott et al. (2011a). In this high-mass sample, slowly rotat- 
ing stars cluster around e(N) = 7.5, which is not explained by 
evolutionary models of the main-sequence. Hence, here again, 
an unidentified mechanism is responsible for the chemical en- 
richment of those objects. We warn about overinterpretating 
the present results. A rigorous comparison of measured nitro- 
gen abundances with theoretical prediction in a Hunter-like dia- 
gram requires population synthesis as performed in Brott et al. 
(201 lb), which we do not perform in this paper. 

Whatever the exact predictions of the models, another con- 
clusion of our study is that on average, the higher the mass, the 
stronger the N-enrichment. This is clear from Fig. 8, in which 
we see that above 100km s , stars with M<27 M have e(A0 
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Fig. 8. Nitrogen abundance (in units of 12+log(N/H)) versus projected 
rotational velocity for the O stars analyzed here. The upper (lower) 
panel shows star swith a mass higher (lower) than 27M . Evolutionary 
tracks for different masses are shown in both panels. They were taken 
from Brott et al. (2011a) and have initial rotational velocities ~100, 
~200, ~300 and ~400 tans'. A horizontal shift corresponding to a 
factor n/4 to account for an average projection effect has been applied 
to the theoretical curves. 



below 7.6, while stars with M>27 M Q have e(N) above 7.9. This 
is a qualitative confirmation of the role of rotational mixing. 



5.3.3. Carbon-to-nitrogen ratio 

In massive stars, the total mass fraction of CNO remains roughly 
constant during the CNO cycle. In the very early phases of evo- 
lution, the CN-equilibrium is achieved first, and nitrogen is en- 
hanced at the expense of carbon.The oxygen content remains 
roughly constant, and somewhat later enters the burning cycle. 
The N/C ratio is an excellent tracer of the chemical evolution 
of massive stars, because it increases substantially when CNO 
cycle-processed material appears at the stellar surface. 

In Fig. 9, we present a comparison between the measured 
N/C ratios (by number) as a function of stellar luminosity. Stel- 
lar evolution models for M„ - 20, 25, 40, and 60 M Q by Brott 
et al. (2011a) are also shown. For each mass, we selected three 
initial rotational velocities (100, 300, and 400 kms -1 ) to cover 
the range of observed present v sin i of our sample stars. The am- 
plitude of the N/C ratio is sensitive to differences in rotational 
velocities, masses, and ages. We can therefore use this plot to 
check whether a star is evolving according to the prediction of 
standard evolutionary tracks. We note that initial abundances 
play a critical role, because the relative increase of the nitrogen 
abundance depends on the initial amount of carbon available for 
conversion into nitrogen, as measured by N/C ratios. We note 
in addition that although the absolute value of the initial C, N, 
and O abundances remains a vigorously debated topic, especially 
in the solar case, the nitrogen-to-carbon abundance ratio (N/C = 



0.13 in the SMC, see Table 3 and Fig. 9) is fairly well established 
and may be considered a better reference point for comparison. 
The direct comparison between the derived N/C and lumi- 
nosities and the predictions of stellar evolution leads to the fol- 
lowing conclusions: 



1. 



M, < 23 M G : four out of the five filled circles, corresponding 
to usin i < 200 kms" , can be accounted for by the 25 M 
evolutionary tracks with rotational velocities between 100 
and 300 km s -1 . For these objects, normal stellar evolution is 
therefore able to explain the observed properties. One slowly 
rotating object (MPG 523) has N/C=5, much higher than pre- 
dicted. A high initial rotation (of about 400 kms -1 ) would 
be needed to reproduce its position in the N/C vs L diagram. 
Of the three fast rotators (open symbols), two are well repro- 
duced by the 400 kms -1 track. The latest object (AzV 267) 
displays a lower enrichment than expected for its luminos- 
ity and v sin ;'. In conclusion, six out eight stars in this mass 
range can be accounted for by stellar evolution. 
23 < M, < 33 M : of the eight objects with usin i < 200 
kms -1 , six have N/C < 1.0, consistent with the 25 M Q evo- 
lutionary tracks with initial rotational velocities below 300 
km s . The two remaining objects (MPG 368 and AzV 446) 
have N/C~2 and log L/L G ~ 5.3. The 25 M evolutionary 
track with usin i = 300 kms -1 can marginally account for 
these stars, provided they are very close to the end of the 
main sequence. Three stars have high rotation rates. Two are 
reproduced by evolutionary models with u ro , - 400 kms -1 , 
one AzV 117 has a luminosity and N/C too high for the 25 
M G track. But its mass (3 1 .9 M ) is intermediate between the 
25 and 40 M track, the latter being able to reproduce its the 
position. Consequently, AzV 117 is marginally reproduced 
by stellar evolution. In total, in the 23-33 M Q mass range, 
eight stars follow normal evolution, and three are marginally 
explained. 

33 < M* < 50 M G : the three objects in this mass range are no 
fast rotators. Two objects (MPG 324 and AzV 243) fall on 
the 40 M evolutionary track with u rot = 300 km s -1 , near the 
end of the main sequence. MPG 324 is a suspected binary, 
which obviously might affect its evolution and surface abun- 
dances. One object (AzV 388) has a higher N/C that cannot 
be explained by moderate rotation. Its usin i is 150 kms -1 , 
higher than that of MPG 324 which has a similar luminos- 
ity, gravity, and mass, but a lower u sin i (70 km s -1 ). Hence, 
AzV 388 might have been initially a fast rotator and has now 
braked to a moderate usin i. In conclusion, the three objects 
in this mass range are marginally explained by normal evo- 
lution. 

M, > 50 M : the only star in this range (MPG 355) is 
marginally explained by the 60 M track at 300 km s -1 if the 
star is close to the end of the main sequence. Given its spec- 
tral type (ON2III(f*)), this is a likely possibility. The nitro- 
gen abundance derived for this star agrees very well with the 
value measured by Rivero Gonzalez et al. (2012) from then- 
analysis of the optical nitrogen lines. These authors found 
values in the range [N]=7.98-8.10 (depending on the tem- 
perature solution), which is consistent with their conclusion 
that their method and tool (the FASTWIND code) for ana- 
lyzing the earliest O-stars would yield slightly lower abun- 
dances than those based on CMFGEN. Nevertheless, the ni- 
trogen abundances from both studies agree within the error 
bars. The N/C ratio is thus much higher than predicted by 
stellar evolution models with moderate initial rotation. Only 
models with v ro , as high as 300 kms -1 can account for the 
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Fig. 9. Nitrogen-to-carbon ratio as a function of luminosity. Circles 
correspond to stars with M < 23 M Q , squares indicate stars with 23 < 
M, < 33 M Q , triangles show stars with 33 < M, < 50 M and pen- 
tagons show stars with M, > 50 M . Open (filled) symbols correspond 
to stars with usin i > 200 kms~' (<, resp.). Stellar evolution models 
taken from Brott et al. (201 la) for M, = 20, 25, 40 and 60 M e are over- 
plotted: dashed (solid, dotted) lines correspond to models with v sin i > 
400(300, 100) km s _1 . 



observed ratio. At SMC metallicity, a 60 M (or above) star 
with such a high rotational velocity may have chemically ho- 
mogeneous evolution (cf. Brott et al. 201 la). 

In total, 15 out of 23 stars (65%) are correctly explained by 
the Brott et al. (201 la) evolutionary tracks; 6 out of 23 (26%) are 
marginally explained, and only 2 objects (9%) are clear outliers. 
The bulk of the SMC O dwarf stars follows normal evolution. A 
fraction of 10 to about 30% of the SMC O dwarfs deviates mildly 
or strongly from normal evolution. For comparison, Brott et al. 
(201 lb) reported that about 40% of the SMC B stars analyzed 
by Hunter et al. (2008) cannot be reproduced by their evolution- 
ary models (the same we used in the present study). Using the 
same dataset, Maeder et al. (2009) concluded that the fraction of 
stars escaping normal evolution was about 20%. Given the un- 
certainty concerning both the value of B stars following unusual 
evolution and the number of true outliers in our sample, we can 
conclude that the majority of SMC O dwarfs are well explained 
by single-star evolution with rotation and that a fraction simi- 
lar to that observed for B stars might need additional physical 
explanation (such as binarity or magnetism). 

5.3.4. CNO 

Achieving the complete CNO-equilibrium requires that hydro- 
gen has been burnt in significant amount in the stellar core. 
Changes in the oxygen surface abundance should therefore be 
observed preferentially in the later evolutionary stages or in the 
more massive and faster rotating stars. In Fig. 10, we com- 
pare measured N/C versus N/O ratios with those predicted by 
stellar evolution models for different masses and rotation rates 



(Brott et al. 201 la). This plot indicates whether the variations of 
N/C and N/O follow the trend expected for material processed 
through the CNO-cycle. The amplitudes of the departure from 
the initial abundance ratios are notoriously sensitive to differ- 
ences in rotational velocities, masses, and ages. They might also 
reflect departures of the initial mixture from the standard CNO 
abundances measured in the SMC interstellar medium (Kurt & 
Dufour 1998), and used in the stellar evolution models. Our re- 
sults can be summarized as follows: 

1. M* < 23 M Q : four out of the five filled circles, corresponding 
to v sin i < 200 km s -1 , can be accounted for by evolutionary 
tracks with rotational velocities between 100 and 300 km s -1 
and masses 20 M and/or 25 M . Their surface properties 
are therefore compatible with the mixing processes as im- 
plement in the evolution models. Again, MPG 523 does not 
follow this trend. The N/C and N/O ratios we measured are 
compatible with a 20 M model with 300 kms~' (but N/O 
would need to be twice as high to be consistent with the 20 
M model and v ro , « 400 kms -1 that was found to best ac- 
counts for MPG 523 behavior, see Sect. 5.3.3). The location 
of MPG 523 in Figs. 9 and 10 strengthens the case that this 
star is indeed fast-rotating. This also implies an evolutionary 
stage rather later than expected from the H-R diagram in Fig. 
4, but consistent with a location on isochrone 5 Myr (see Fig. 
5). 

Of the three fast rotators (open symbols), two are well repro- 
duced by the 400 km s" 1 track, these are AzV 267 and MPG 
299. The other object (NGC 346-046) displays an oxygen 
depletion much too weak for the measured carbon and nitro- 
gen abundances. To match the prediction of a 20 M and 
v rot ~ 400 km s -1 model, O/H would need to be roughly four 
times lower than measured. Although the number of oxygen 
lines for abundance determination is limited and the sensi- 
tivity/accuracy reduced because of the high v sin i of the star, 
we can rule out such a low abundance. At this point, six out 
eight stars in this mass range can be accounted for by stellar 
evolution. 

2. 23 < M* < 33 M : Seven out of eight objects with usin i 
< 200 km s _I display N/O and N/C ratios consistent with the 
25 M evolutionary tracks with initial rotational velocities of 
about 100 kms -1 . For the one remaining object (MPG 368), 
this evolutionary track is able to account marginally for its 
surface abundance but implies that the star is at the terminal 
stage of its evolution (the end of helium-core-burning in this 
grid), which can be ruled out from its other observed proper- 
ties. 25 M models with faster rotation (u sin i > 300 km s" 1 ) 
can better account for the surface abundances of MPG 368. 
This behavior is easier to reconcile with the likely binary sta- 
tus of this object. 

Of the three stars with high rotation rates, two are reproduced 
by evolutionary models with 25 M and i>sin i = 300-400 
kms" 1 . The last one (AzV 177) has an N/C too high for 
any 25 M Q track, but N/C and N/O are compatible with a 
model with 30 M and v rot « 300 kms -1 , consistent with 
the conclusion reached in Sect. 5.3.3. Including the present 
result, we argue that AzV 1 17 is indeed reproduced by stellar 
evolution. In total, in the 23-33 M mass range, all (but one) 
stars follow normal evolution. 

3. 33 < M t < 50 M : Two out of three stars in this mass range, 
namely AzV 388 and AzV 243, fall on the 40 M evolution- 
ary track with v ro , w 400 kms -1 (AzV 243 is also on the 40 
M , 300 kms -1 track), although none of them is a fast ro- 
tator. MPG 324, on the other hand, is a clear outlier in this 
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diagram for this mass range, possibly related to its binary 
status. MPG 324 and AzV 388 have very different N/O and 
N/C ratios, although they have similar luminosity, gravity, 
and mass, but their v sin i differs by a factor of two. Together 
with the behavior observed in Fig. 9, this reinforces the con- 
clusion that AzV 388 initially was a fast rotator and has now 
braked down to a moderate v sin i. 
4. M, > 50 M Q : As the only star clearly more massive than 
50 M G (but whose real mass is unknown), MPG 355 has 
N/O and N/C ratios consistent with the 60 M G tracks hav- 
ing u rot « 300 kms -1 . This again strengthens our interpre- 
tation throughout this paper that this star is a very massive 
fast-rotating object in a chemically homogeneous evolution. 

In summary, 19 out of 23 stars have surface abundance pat- 
terns that are compatible with normal stellar evolution as imple- 
mented in the Brott et al. (201 la) evolutionary tracks. Only four 
objects are clear outliers, one of them is marginally consistent 
with the expected evolution for the correct mass range although 
of unlikely evolutionary status. 

For the stars that deviate from stellar evolution predictions, 
the amplitude of the deviations are generally large enough that 
differences in the initial abundances could not help to reconcile 
the outputs from stellar evolution models and our measurements. 
It is indeed well-known that the slope of the theoretical tracks 
is very sensitive to the adopted initial mixture (see for instance 
Fig. 3 of Przybilla et al. (2010) and a span of N/O and N/C ra- 
tios is expected for different initial chemical composition. In our 
case, all outliers but one (MPG 523) have C+N+O in good con- 
sistency with the values measured in the H II regions (Kurt & 
Dufour 1998). We furthermore stress that several stars that fol- 
low the normal evolution do present some departures from these 
ISM values, on the other hand . However, the scatter around the 
SMC abundances is small for most objects, which indicates that 
the initial composition of the stars is compatible with the ISM 
values. This rules out a difference from standard ISM values in 
the initial abundances as the origin for the observed behavior of 
the outliers stars. 

Alternatively, including magnetic fields in stellar evolution 
models could change the speed and amplitude with which abun- 
dance patterns are built up at the surface. A recent study by 
Meynet et al. (201 1) confirmed that including magnetic braking 
yields very different results for the surface abundances, which 
depend on the assumed rotation law inside the star. For stars 
with differential rotation, mixing is faster and stronger, but for 
stars with solid-body rotation, mixing is inhibited because the 
stars spin down rapidly. In this case, surface abundance ratios are 
lower than in models without magnetic braking. Alternatively, 
(close) binarity would affect the evolution and likely the surface 
abundances of the stars (from mixing and possibly mass trans- 
fer), although it would also imply different luminosities, hence 
evolutionary masses. However, a comparison with models tai- 
lored for magnetic braking (and/or lower masses) would be re- 
quired for more quantitative conclusions. 

The overall properties of the sample stars, including their lo- 
cation in the H-R diagram and surface abundances, indicate that 
they are observed at different stages of their chemical evolution, 
even though they are of very similar age. The underlying drivers 
of the observed differences are unclear and could be manyfold - 
from enhanced mixing related to binary evolution, or to evolu- 
tion with rotation and a magnetic field (Meynet et al. 201 1). 

Recent results reported in Sana et al. (2012a,b) suggest that 
the binary fraction among massive stars is at least as high as 
50%. Most of the (potential) binaries we identified follow a nor- 
mal, single-star-like, evolution. These are of the late spectral 



type and are expected to evolve more slowly, which might ex- 
plain that no signs of binary evolution is visible yet. On the 
other hand, we note that the most two massive binaries, namely 
MPG 324 and MPG 368, do present surface abundance proper- 
ties at odds with normal stellar evolution. Two circumstances 
may boost this behavior. First, because they are more massive, 
these stars are expected to evolve faster and display signs of 
chemical evolution more quickly. Second, they are potentially 
close-binary systems (we recall that radial-velocity variations 
have been found for these stars); this is very likely to impact 
the evolution of these objects even more. 



5.4. Masses 

Spectroscopic masses have been obtained from line profile fit- 
ting, while evolutionary masses have been derived by bilinear 
interpolation on the tracks and isochrones in the H-R diagram. 
Since model atmospheres only account for radiation pressure 
plus gas pressure to balance gravity, spectroscopic masses are 
correct only if the centrifugal acceleration is negligible, that is, 
only for slow rotators. Otherwise, spectroscopic analyses are 
expected to underestimate the stellar mass. For this reason, we 
corrected the spectroscopic surface gravities for centrifugal force 
induced by rotation, following the approach outlined in Repo- 
lust et al. (2004). Assuming a random distribution of orienta- 
tions for the rotational axis, the corrected surface gravity is the 
sum of the effective gravity and a correction term where the cen- 
trifugal acceleration averaged over the stellar disk is given to a 
very good approximation by the projected centrifugal velocity: 
g c - g eff + (v sin i) 2 /R t , 

where g c and g e fj are the centrifugally corrected and effective 
surface gravity, respectively. In most cases, the correction is 
small, with A log g < 0.05 dex, implying corrections of 10% 
or less of the spectroscopic masses. Only stars with usin i > 
300 km s- 1 - NGC 346-046 and MPG 299 - are affected by a 
larger correction of 0.06 dex and 0. 1 dex respectively, implying 
a correction of the stellar mass as high as 25% in the later case. 

A comparison of spectroscopic and evolutionary masses as 
listed in Table 2 is displayed in Fig. 11. The errors for the evolu- 
tionary mass reflect the errors spanned by the stellar luminosity 
and effective temperature, respectively, for the adopted set of 
evolutionary tracks, i.e., for a given initial rotation rate. With 
the introduction of rotation in stellar evolution models, it is no 
longer possible to assign a unique value for the stellar mass from 
a given set (L*, T e g), and therefore this question translates into 
larger uncertainties on the evolutionary masses. Likewise, the 
large error bars on spectroscopic masses result from uncertain- 
ties on log g. Since the stellar distance is known within a fairly 
good accuracy for SMC stars, the error in \ogg is indeed the 
critical source of error in the spectroscopic mass, hence in the 
ratio of the evolutionary mass to the spectroscopic mass. The un- 
certainties in the determination of surface gravities are partially 
rooted in the use of echelle spectra, whose rectification is diffi- 
cult, especially in the vicinity of broad H lines. The accuracy 
of the surface gravity measurements is considerably improved, 
however, when optical spectra can be used to constrain log g in 
the models (see also Sect. 3.2). This is the case for the 12 stars 
in NGC 346 and for AzV 177, AzV 388 and AzV 243. Four 
more stars have 2dF spectra, that yield a relatively accurate de- 
termination of log g despite a lower spectral resolution. For the 
five stars with UV spectroscopy only, we adopted a higher error 
bars on log g that is indicative of log g scatter for dwarf O stars 
around the assumed value of log g- 4.0. 
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Fig. 10. N/C ratio as a function of N/O. The same symbols and color coding as in Fig. 9 have been adopted. Stellar evolution models according to 
Brott et al. (201 la) for M, = 20, 25, 40 and 60 M Q are overplotted. Here again, the same line styles as in Fig. 9 were adopted; we also color-scoded 
each corresponding mass range for clarity. The upper left panel presents the complete sample, the three remaining panels display the three mass 
ranges discussed in the text. 



Figure 1 1 shows that, within the error bars, most stars follow 
the 1:1 relation, although three objects still present (i.e. even 
within the error bars) the mass-discrepancy first pointed out by 
Herrero et al. (1992). A fourth object is also only marginally on 
the 1:1 relation. The study by Mokiem et al. (2006) also con- 
cluded that a mass-discrepancy was present in some SMC stars. 
They argued that the discrepancy observed in helium rich stars of 
their sample was possibly a consequence of fast rotation in these 
objects. Our work shows that one star showing the discrepancy 
is indeed helium rich, namely MPG 487 but MPG 299, who is 
the other helium rich star and is a fast rotator, shows a very good 
agreement between the spectroscopic and evolutionary masses. 
Furthermore, the other stars showing the mass discrepancy are 
easily identified as suspected binaries (AzV 326, AzV 189, AzV 
148) and are not helium enriched objects (hence not likely to 
rotate rapidly). We note that for MPG 355 whose evolutionary 
masse is not constrained from the grid by Brott et al. (2011a), 
the mass discrepancy could be large, as pointed out in Rivero 
Gonzalez et al. (2012). 

A different view of our results is provided by the histogram 
of the mass ratios M spec / M evo i, which we present in Fig. 12. 
The distribution of the ratio divides into two peaks. In total, five 
stars present mass ratios higher than 1.2. Four of these stars are 
clear binaries (AzV 326, AzV 189, AzV 148, and MPG 487) and 
their masses, hence M spec / M evo i ratios, are affected by larger un- 



certainties. We also note that for MPG 368, also identified as a 
binary (see also Massey et al. 2009), the mass ratio deviates from 
the 1 : 1 relation by the same amplitude as the four later-type bina- 
ries. The bias toward highly discrepant ratios for these five stars 
is therefore probably related to the binary status of these objects. 
We note, on the other hand, that the earliest-type binary, MPG 
324 instead presents an M spec / M evo i ratio that excellently agrees 
between evolutionary and spectroscopic masses. Furthermore, 
we were not able to identify any issue in the spectroscopic de- 
termination of the surface gravity of the fifth star, namely MPG 
682. A lower log g than quoted in Table 2 is firmly excluded; 
this would decrease the spectroscopic mass and help reconcile it 
with the evolutionary mass. 

For the large majority of the sample stars, the distribution is 
skewed toward values lower than unity (with a mean of 0.92 and 
a standard deviation of 0.08). This indicates that even includ- 
ing the corrections for rotation applied to spectroscopic masses, 
the latter are overall smaller than evolutionary masses. Obvi- 
ously, this comment holds for the adopted set of evolutionary 
tracks only. Evolutionary masses (and ages) would be different 
if we adopted tracks and isochrones with different initial rota- 
tion rates. The lack of knowledge of the actual rotation rates of 
the stars hampers the use of evolutionary models computed with 
rotation rates appropriate for each individual object.This trans- 
lates into larger uncertainties on the actual evolutionary masses 
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Fig. 11. Spectroscopic versus evolutionary masses. The red filled 
squares show the rotators with usin i > 300 kms~' (see text). Green 
symbols indicate masses as determined in Mokiem et al. (2006) for the 
stars in common with the present study. MPG 355 is not plotted because 
its evolutionary mass could not be determined with the model grid we 
used (Brott et al. 2011a). 

and the observed skewness of the histogram in Fig. 12 probably 
reflects this effect. Quantitatively, however, we found that M evo ] 
for our target stars would be changed by half a solar mass on 
average if we used evolutionary models with v rot ~ 100 kms -1 
for instance. With these models, though, the youngest two stars 
(NGC 346-046 and MPG 523) of our sample would be located to 
the left of the ZAMS, even within the error bars. For such cases, 
estimates of evolutionary masses are highly uncertain. We pre- 
sented arguments in Sect. 5.2.1 that showed that NGC 346-046 
is potentially best described by a model for a 20 M Q star rotat- 
ing between 400 km s -1 and 500 km s (also consistent with the 
chemical N/C pattern, see Sect. 5.3). The derived evolutionary 
mass would be roughly 3 M G lower than quoted in Table 2, a per- 
fect match with the spectroscopic mass. We stress furthermore 
that the same models also account very well for the position of 
MPG 523 in the H-R diagram, while providing a better match to 
the measured N/C ratio. 

6. Summary 

We have modeled the UV and optical spectra of 23 O-type dwarf 
stars in the SMC to derive their surface and evolutionary proper- 
ties. Our results can be summarized as follows: 

1. We can derive accurate luminosities from modeling the 
SEDs of the sample stars. We found three objects with lu- 
minosities too high for their spectral types. We showed that 
these stars are very likely binaries consisting of two objects 
of the same or very similar spectral types. Three other stars 
were already identified as potential binaries; we confirmed 
one of them from SED fitting. The measured fraction of bi- 
naries of our sample is 26%, which agrees well with more 
systematic dedicated studies of binarity in massive stars. 

2. The sample consists of stars in the giant H II region NGC 346 
and field stars. From the H-R diagram, we found that the 
younger stars belong to NGC 346, while fields stars are 
2 Myr old or older. The stars in NGC 346 have ages rang- 
ing from younger than 1 Myr to about 6 Myr (for MPG 012, 
the only subgiant of this sample). 



Fig. 12. Histogram of the mass-ratio distribution. A bin size of 0.1 is 
adopted for the plot. Here again, MPG 355 is not included (see text for 
comment). 



3. We found two fast rotators with osin i > 300 kms ', 
MPG 299 and NGC 346-046. Both stars are 1 Myr old at 
most, as indicated by their location on or very close to the 
ZAMS. The offset in the position of these two fast rotators 
compared to the other stars confirms the predictions of evo- 
lutionary models that fast-rotating stars tend to evolve more 
vertically in the H-R diagram. 

4. Three stars of our sample, MPG 113, MPG 523 and 
NGC 346-031, were originally classified as luminosity-class 
Vz. This luminosity class points to massive stars whose 
spectral properties are expected to indicate extreme youth. In 
contrast to Mokiem et al. (2006), our analysis indicates that 
only MPG 523 is on the ZAMS, the other two stars are more 
evolved, even within the error bars. For MPG 113, a dif- 
ference in the effective temperature of 3,000 K between both 
studies is the cause of the shift in age. The case of NGC 346- 
03 1 is very interesting because it points to a different inter- 
pretation of the strength of the He n /14686 line, the defining 
criterion of the Vz class, based on different mass-loss rates 
and helium abundances. 

5. The discrepancy between spectroscopic and evolutionary 
masses is reduced but still present. For a subsample of 17 
of the 23 stars, the mean mass ratio is 0.92, with a standard 
deviation of ±0.08. Binary stars present stronger deviations 
from the 1 : 1 relation of M spec / M evo i (i.e., a deviation of 20% 
or more), while one single star is the only clear outlier to 
the M spec = M evo ] relation. Obviously, the detailed numeri- 
cal results of this analysis are valid only for the adopted set 
of evolutionary tracks and isochrones. Investigation of the 
impact of different rotation rates on the evolutionary calcu- 
lations showed rather limited influence on the mass determi- 
nation, however, not enough to explain the observed mass 
discrepancy. 

6. Photospheric UV lines of C and N were used to derive sur- 
face abundances. Our results may be regarded as an ex- 
tension toward higher masses of the work by Hunter et al. 
(2009) on B-type stars. Using their published abundances 
for B stars, we found that the distribution of O and B stars in 
the e(N) - v sin i diagram is the same, especially for the low- 
velocity part of this diagram. This suggests that the mecha- 
nisms responsible for the chemical enrichment of slowly ro- 
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tating massive stars depends weakly on the mass and that O 
stars are at least as enriched as B stars. 

7. We compared the measured N/C ratios as a function of stellar 
luminosity, and the predictions of the stellar evolution mod- 
els by Brott et al. (2011a). Sixty-five percent of our sam- 
ple stars are well reproduced by evolutionary models includ- 
ing rotation. Twenty-six percent are marginally compatible 
with these models, and 9% are clear outliers. The fraction 
of O dwarfs in the SMC that cannot be explaned by standard 
evolution is similar to that of B stars, considering the uncer- 
tainties in the estimate of the exact fraction of these outliers 
(between 20 and 40%). 

8. Finally, a majority of stars, irrespective of their mass, have 
N/C ratios as a function of N/O that match the predicted be- 
havior well, which strengthens the physical assumptions for 
mixing that are currently implemented in the evolution mod- 
els presented in Brott et al. (201 la). Two out of the four stars 
that differ from this normal evolution are massive, possibly 
close, binaries, and their surface abundances, should be in- 
vestigated with ad hoc evolutionary models. The fraction of 
stars close to the ZAMS that show unexplained differences to 
standard evolution is thus 9%. This fraction is significantly 
smaller than the a; 40 % outliers of B-type core-hydrogen- 
burning stars in the SMC found by Hunter et al. (2009) and 
Brott et al. (2011b). Whether or not this reflects the actual 
trend for O stars needs to be investigated by extending the 
present analysis to more evolved objects in the giant and su- 
pergiant luminosity classes, such that the compared samples 
of O- and B-type stars are as similar as possible. 

The different conclusions that we reached regarding the ni- 
trogen surface abundance on one hand and the N/C or/and N/O 
ratios on the other hand deserves more investigation. Maeder 
et al. (2009) and Brott et al. (201 lb) previously pointed out that 
the so-called Hunter diagram might not be the best way to con- 
strain rotational mixing. These authors discussed how the rela- 
tive influences of binarity, mass loss, metallicity, age, and ro- 
tation on the chemical evolution of massive stars could blurr 
the interpretation of this diagram. Our results also support this 
conclusion. We found indeed that a (very) limited fraction of 
stars do not comply with the predictions of evolution models 
with rotation when abundance ratios are used. Taken at face 
value, they suggest that the elements abundances taken sepa- 
rately might not be as relevant as abundance ratios when obser- 
vations are to be compared with outcomes from stellar evolution 
models. In other words, the more constraints are obtained on 
surface abundances (beyond N abundance only), the better the 
agreement between models and observations. A larger sample 
including more evolved objects would allow one to investigate if 
these results vary along the evolution. We will address this ques- 
tion in a forthcoming paper where a similar spectral analysis of 
SMC giant and supergiant stars will be presented. 
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Appendix A: Atomic data 

We provide below a brief description of the atomic data used in 
the present set of calculations. For many atomic species differ- 
ent atomic data sets are available. In the present calculations we 
did not always use the latest data set, although test calculations 
show that the choice of other data sets would not affect our con- 
clusions. 

For H we used collision strengths adapted from Mihalas et al. 
(1975), which are similar to those used in tlusty (Hubeny & 
Lanz 1995b). For the Stark broadening of optical hydrogen 
lines we used the results from Lemke (1997); for the optical 
lines they give similar results as are found in Schoning & But- 
ler (1989c). Collision rates for Hei were taken from Berring- 
ton & Kingston (1987), while photoionization cross-sections 
are "below-resonance" fits to the calculations of Fernley et al. 
(1987). For Hen we used the Stark- broadening calculations of 
Schoning & Butler (1989c,b,a). For He 1 ,14471 and i4922 we 
used data from Barnard et al. (1974), for i4389 and i4026 the 
data was taken from Shamey (1969) (kindly supplied by Keith 
Butler), while for other He i lines we used data from Dimitrijevic 
& Sahal-Brechot (1984). Other data for H and He are discussed 
in Hillier (1987). 

Civ oscillator strengths were taken from Leibowitz (1972) 
while the photoionization data were taken from Peach et al. 
(1988) and Leibowitz (1972). C m oscillator strengths, photoion- 
ization data, and dielectronic data were taken from P. J. Storey 
(private communication), while collision rates for the lowest six 
terms were taken from Berrington et al. (1985). 

Nv oscillator strengths for n < 6 were taken from Liang 
& Badnell (2011), the photoionization data were taken from 
Peach et al. (1988), and collisional data for the 2s-2p, 2s-3 
and 2p-3 transitions were taken from Cochrane & McWhirter 
(1983). Niv oscillator strengths and photoionization cross- 
sections were computed in Tully et al. (1990) as part of the 
opacity project (Seaton 1987; Cunto et al. 1993). The identifica- 
tion of the 2s 4s 3 S and 2p 3p 3 S states were switched (see Allard 
et al. 1991). Oscillator strengths for Niv forbidden (and semi- 
forbidden) lines were taken from Nussbaumer & Storey (1979b), 
while for oscillator strengths for low-temperature dielectronic 
recombination transitions, we followed Nussbaumer & Storey 
(1983, 1984) (private communication). Transition probabilities 
for the N m intercombination lines were taken from Nussbaumer 
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& Storey (1979a) which agree reasonably well with the experi- 
mental values found in the recent study of Trabert et al. (1999). 

Ov oscillator strengths were taken from (Nussbaumer & 
Storey 1983, 1984, private communication). Collision rates for 
the six lowest terms of Ov were taken from Berrington et al. 
(1985). Oiv oscillator strengths were computed as part of the 
opacity project (Seaton 1987) and were obtained from topbase 
(Cunto et al. 1993). Intercombination data for Oiv is from the 
compilation of Mendoza (1983). Oiv photoionization data were 
taken from Nahar (1998) and were obtained through norad 4 . 
O iv collision rates were taken from Zhang et al. (1994). O m os- 
cillator strengths and photoionization data were taken from Luo 
et al. (1989). Collisional data for the four lowest terms in Oiii 
were taken from the compilation of Mendoza (1983). 

Oscillator strengths for Ne in and Ne iv were taken from the 
compilation of Kurucz & Peytremann (1975) and the opacity 
project (Seaton 1987). The Neni and Neiv photoionization 
cross-sections were obtained from the opacity project (Seaton 
1987; Cunto et al. 1993). Collisional data for the four lowest 
terms in Nev were taken from Lennon & Burke (1994), colli- 
sional data for the three lowest terms of Ne iv were taken from 
the compilation of Mendoza (1983), while collisional data for 
Ne in were taken from McLaughlin & Bell (2000). 

Oscillator strengths and photoionization cross-sections for 
N m, Ne v (Luo & Pradhan 1989), Si in (Butler et al. 1993), Si iv, 

5 m (Nahar & Pradhan 1993), and S iv were obtained from top- 
base (Cunto et al. 1993), and were computed as part of the opac- 
ity project (Seaton 1987). Intercombination data for Siv were 
taken from the compilation of Mendoza (1983) while collision 
strengths between the 2 P" and 4 P levels were taken from the 
compilation of Mendoza (1983). Collisional data for the four 
lowest terms for S in were taken from Tayal (1997). 

Oscillator strengths for Fe iv, Fe v, and Fe vi were computed 
by Kurucz (2009) and were obtained (before 2001) from Bob 
Kurucz's website 5 . The oscillator strengths of two Feiv lines 
that overlap with the He i resonance line near 58.4 nm were re- 
duced by a factor of 10 (see Najarro et al. 2006). Photoionization 
cross-sections for Fe iv, Fe v, Fe vi and Fe vn were computed as 
part of the opacity project (Seaton 1987) and were obtained from 
topbase (Cunto et al. 1993). Feiv collisional data is from Zhang 

6 Pradhan (1997). 

Appendix B: Best fits 

In this appendix we give more details of the parameter determi- 
nation of each individual star and compare our best-fit models to 
the optical and UV spectra. Absolute fluxes in UV spectra are 
provided after correcting for reddening (see Sect. 3.2 and Table 
1) and are expressed in ergs . A -1 . s _1 . cm -2 

Appendix B.1: MPG 355 - ON2 III (f ) 

MPG 355 was already studied extensively in the UV and optical 
in Bouret et al. (2003), to which we refer for a thorough discus- 
sion of the STIS spectrum. More recently, Massey et al. (2009) 
obtained a new optical spectrum with a higher S/N of MPG 355; 
they detected a very weak He AAAll (their Fig. 7) that is not 
visible in our own spectrum. This biased Bouret et al. (2003) 
toward a higher T e s . 

On the other hand, He i /15876 is not detected in our high 
S/N UVES spectrum while this line would be expected at r e ff 



= 50,000 K in a test model we computed. Furthermore, several 
lines of CNO elements argue in favor of a higher effective tem- 
perature T e ff — 52,500 K. This is the case of C iv Al 169 and C in 
A 1176, where the intensity of individual lines (the C m line is 
not detected in the STIS spectrum) is best fitted for such a T e g. 
Other lines include N v X14604-4620 and N iv ,16381, where 
both the intensity of individual lines and line ratio of these suc- 
cessives ionization stages argue for the adopted r e ff. Finally, 
we found that emission from incoherent electron scattering in 
the wings of He n /14686 is predicted somewhat too strong and 
that a higher effective temperature (55,000 K) would solve this 
problem. Overall, we are quite confident that the effective tem- 
perature of MPG 355 as quoted in Bouret et al. (2003) and in the 
present paper is correct. The recent results by Rivero Gonzalez 
et al. (2012) also strengthen this conclusion and the other quan- 
tities derived here. 

Compared with our previous study, we only changed the ni- 
trogen surface abundance to better fit the N iv /16381, still main- 
taining a very good fit to other lines from nitrogen ions. 



Appendix B.2: AzV 1 77 - 04 V 

Detailed modeling of the optical spectrum of AzV 177 was pre- 
sented by Massey et al. (2005) and Rivero Gonzalez et al. (2012). 
Significant differences are found from these studies that concern 
photospheric quantities such as the surface gravity or the helium 
abundance, and the wind parameters. The 0.2 dex difference in 
log g c is very surprising since we used the same optical spectrum 
(kindly provide by P. Massey) and derived a similar r e ff. Our 
modeling of this star relies on the FUV COS spectrum, which 
offers a better grasp on the wind properties. We also have several 
CNO lines at our disposal to constrain the surface abundances. 
The line ratio C iv Al 169 to C in Al 176 argues for a T^ slightly 
lower than in Massey et al. (2005) and Rivero Gonzalez et al. 
(2012), although both estimates would agree within error bars. 
The mass-loss rates also match within a factor of two, which 
is a very good agreement given the clumping filling factor we 
adopted. The wind velocity exponent is slightly higher in our 
case (J3 - l.l rather than 0.8), a result mostly driven by the fit to 
the C iv resonance doublet. We note that this higher ft addition- 
ally precludes using He/H=0.15 from Massey et al. (2005) and 
Rivero Gonzalez et al. (2012), because He n 4686 would show 
strong emission wings (already too strong here), which are not 
observed. 



Appendix B.3: AzV 388 - 04 V 

AzV 388 was included in the sample of SMC stars investigated 
in Mokiem et al. (2006). We used the same optical spectrum, 
provided by P. Crowther. The fit to this spectrum is excellent 
throughout the whole spectral range, including UV. The derived 
parameters agree very well with those listed in Mokiem et al. 
(2006), well within the error bars. The only major discrepancy 
is the mass-loss rate, which we found to be a factor 3.3 times 
lower than their value which is based mostly on Ha. Still, the 
derived clumping filling factor / = 0.11 implies that our mass- 
loss rate would be scaled up by the exact same factor, if the wind 

4 http://www.astronomy.ohio-state.edu/~nahar/nahar_radiativeatomicdatayVere assumed to be homogeneous. This last assumption is in- 

5 http://kurucz.harvard.edu compatible with the O v /11371 line profile, however. 
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Fig. B.l. Left: Best-fit model for MPG 355 (red line) compared with the COS spectrum (black line). Right: Best-fit model for MPG 355 (red 
line) compared with the AAT and UVES spectra (black line). 
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Fig. B.2. Left: Best-fit model for AzV 177 (red line) compared with the COS spectrum (black line). Right: Best-fit model for AzV 177 (red line) 
compared with the CTIO spectrum (black line). 
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Appendix B.4: MPG 324 - 04 V ((f)) 

MPG 324 was studied in the UV and optical (AAT spectrum) in 
Bouret et al. (2003) and Heap et al. (2006) and was more re- 
cently included in two other studies by Mokiem et al. (2006) and 
Massey et al. (2009), both of them using the code FASTWIND 
for the modeling of the optical spectrum. Using this time the 
same optical spectrum (FLAMES) as in Mokiem et al. (2006), 
supplemented with our own UVES spectrum, the fundamental 
and wind parameters we derived differ only moderately from our 
initial study. r e ff was increased by 600 K (a change driven by the 
He 1 5876 line observed with very good S/N ratio with UVES and 
no apparent nebular contamination) and agrees reasonably with 
other estimates for this star. The only significant difference with 
published results is the projected rotation velocity, which we find 



to be two times lower (usin i = 70 kms _I vs 150 kms~') than 
quoted Massey et al. (2009). Such a high value seems incon- 
sistent with the UV spectrum, however, especially with the C in 
A 1176 multiplet. We finally note that from small velocity shifts 
observed in some lines, Evans et al. (2006) concluded that MPG 
324 was probably a single-lined binary. 



Appendix B.5: MPG 368 - 05 V ((f+)) 

MPG 368 was studied in the UV and optical (AAT spectrum) in 
Bouret et al. (2003) and Heap et al. (2006). Based on their own 
optical spectrum, Massey et al. (2009) confirmed the interpreta- 
tion of Heap et al. (2006) that this star is a binary. A new optical 
spectrum obtained with UVES shows the He i lines as compos- 
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Fig. B.4. Left: Best-fit model for MPG 324 (red line) compared with the STIS spectrum (black line). Right: Best-fit model for MPG 324 (red 
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ites of two components (see Fig. B.5). Higher ionization lines 
seem to be unaffected, either in UV or optical, suggesting that 
the other component of this system is a late-O or B-type star. 
The fit to the observations with our model proves to be very rea- 
sonable, confirming that the companion of MPG 368 contributes 
only marginally to the wind line spectrum. 

Appendix B.6: MPG 113- OC6 Vz 

MPG 113 was extensively studied in Bouret et al. (2003), Heap 
et al. (2006), and Mokiem et al. (2006). Our best-fit model has 
the same parameters as quoted in our first study but for the mass- 
loss rate, which we slighlty increased to 3x 10~ 9 M Q yr~' as a 
result of an improved log Lx/L\, \ ratio used in the new model 
(too high in the former study). Our T & g is still 3000 K lower 
than found in Mokiem et al. (2006), and although we used the 
same FLAMES spectrum as they did, we did not need to increase 
the helium abundance to fit the helium lines. This difference 
might come from a bias toward lower temperatures to reproduce 
the observed ratio of lines C iv Al 169 / C in Al 176 in our case. 
For the adopted carbon abundance (based on these UV lines), 
the optical lines of C in AA4647 -50-52 were predicted slightly 
too strong, while C hi /15696 was slightly too weak. Both are 
nevertheless predicted in emission, as observed. We note that 
the N m /L14634-4640 lines are neither observed nor predicted in 
our model. The resulting nitrogen-to-carbon ratio is above solar 
(see Table 3), like other stars of the 06 class, and as discussed 
in Sect. 5.2, the reason for the OC6 classification of MPG 1 13 is 
more complex than a simple abundance criterion. 



Appendix B. 7: AzV 243 - 06 V 

Of the 06 V stars of our sample, AzV 243 has the largest lumi- 
nosity, consistent with its photometry. The luminosity we derive 
agrees well with the one quoted in Mokiem et al. (2006), but we 
found this star to be significantly cooler (by 2500 K) than in this 
study. With its higher luminosity, AzV 243 does not fall into the 
weak-wind regime, which is also confirmed by its COS spec- 
trum, showing clear PCygni profiles at N v and C iv resonance 
doublet. The wind acceleration parameter /3 - 1 is intermediate 
between the value /3 - 0.7 given in Haser et al. (1998) and the 
value /3 - 1.37 in Mokiem et al. (2006). The latter is prone to 
higher uncertainty, however, because it was derived from the Ha 
line mostly, which poorly responds to wind parameters for M < 
1.x 10 7 M G yr~' as is the case here. 



Appendix B.8: AzV 446 - 06.5 V 

AzV 446 was studied in Massey et al. (2005). Based on the line 
widths of He i /14471 and He n /I4542, these authors noted that 
this star could be a binary although the absolute magnitude of 
AzV 446 was consistent with a single star of luminosity class V. 
The very good agreeement of their model with the observed op- 
tical spectrum also supported the single-star interpretation. Our 
modeling of the flux-calibrated UV spectrum of AzV 446 addi- 
tionally confirms this (the SED of our model also matches the 
photometry up to NIR). The fundamental parameters of our best 
model marginally differ from those in Massey et al. (2005) but 
for the wind parameters. From the UV only, though, we were 
unable to constrain the surface gravity and helium abundance to 
a good accuracy and simply adopted standard values for a young 
main-sequence star. We note that Massey et al. (2005) measured 
y = 0.15, surprisingly high for such a star. From the very weak 



C iv resonance profile, M - lx 10~ 7 Mq yr -1 is ruled out. How- 
ever, this profile, although weak, is clearly seen, as is the N v 
AA 1248- 1250 profile, showing a discrete absorption component 
near the blue edge of the line. We also note that the C iv AA 1548- 
1550 resonance line presents a blueward extension that is clearly 
sensitive to /. We find that these two wind profiles are best re- 
produced when M - 4x 10~ 9 M G yr _1 , / = 0.1, and Voo = 1400 
kms~'. 



Appendix B.9: MPG 356 - 06.5 V 

This star was modeled in Heap et al. (2006), although in the 
UV only. Using the UVES spectrum we obtained for this star, 
we derived a higher temperature and surface gravity than Heap 
et al. (2006), although still within their error bars. Our lumi- 
nosity is also lower (0.4 dex) but consistent with the complete 
photometric dataset, from UV to NIR. From the photospheric 
lines of carbon and nitrogen in the UV spectrum, we derived 
surface abundances significantly different from those in Heap 
et al. (2006). The carbon abundance is higher than measured 
in H II region of the SMC or even than a solar carbon abundance 
scaled down by a factor 0.2 corresponding to the global SMC 
metallicity. This abundance is consistent with the UV lines and 
with the C in /15696 line in the UVES spectrum. Heap et al. 
(2006) found a moderate (if any) nitrogen enhancement from 
N in AAl 182-1 184, while we derived an enrichment by up to 
a factor of seven (compared to the adopted baseline), from the 
same STIS spectrum. The N/C ratios, on the other hand, agree 
in both studies. Keeping the carbon and nitrogen abundances 
fixed, we constrained the mass-loss rate from N v AA 1248- 1250 
mostly, because the C iv resonance doublet shows no clear wind 
profile. An upper limit on the mass-loss rate M < 3.5x 10~ 9 
M Q yr~' was found. 

Appendix B. 1 0: AzV 429 - 07 V 

AzV 429 is substantially more luminous than the other 07 V 
stars of the sample (see below and Table 2). It is also probably 
more evolved, as indicated by its position in the H-R diagram. 
Having only the UV spectrum for this star, our estimate of the 
surface gravity is based on the Fe v to Fe iv index from Heap 
et al. (2006) and suffers from higher uncertainty. The surface 
abundances of carbon and nitrogen are such that the N/C ratio is 
one of the lowest of the whole sample and only twice the inter- 
stellar value. Nitrogen is enhanced, but the carbon abundance is 
higher than listed in Kurt & Dufour (1998). The C iv ^1548- 
1550 doublet presents a weak but well-defined wind profile. The 
best fit to this profile is obtained using M - 4.0x 10~ 9 M Q yr~' 
and / = 0.05. 



Appendix B.11: MPG 523 - 07 Vz 

Our UVES spectrum of MPG 523 is not as strongly contami- 
nated with nebular emission as was the FLAMES spectrum used 
in Mokiem et al. (2006). Both spectra can be best fitted with r e jf 
= 38,700 K and a relatively high surface gravity, logg = 4.25. 
The fit to the helium lines is good except for He i /15876, which is 
too broad in our model, while the observed profile also presents 
a warp on its blue side. The usin i is moderate (50 kms~') and 
consistent with the UV spectrum as well. This star is not listed 
as a potential or known binary, nor do we find it overluminous 
for its spectral type. We were unable to find a satisfactory ex- 
planation for this profile. The C iv /LU548-1550 shows up as 



Article number, page 28 of 43 



Bouret et al.: Massive stars at low metallicity 




5.0x10"" - 



8.5x10"' 
2.0x10"' 
1.5x10"' 
LOxlO"' 



00 

S-, 

><; 

3 

3 

o 

w 

< 



Wavelength (E) 



X 



O 

E 




5620 5660 




50 


1 158 




1167 


1175 




1183 






1192 12< 


\ 


1 


1 1 I 


1 


/* 












^WPfTfl 


i 


,./ 




1 


1 1 




1 



DO 


1242 


1263 


325 


1367 


140B 


14 


l\\ 


1 
1 




1 
1 


1 1 







E 
-a 



~ 



j q I »i w M i wyw ■ 



.. * ■*■■■■» .ywiHn*^. 







"> rt >r ■ ^Um %» v -/ w ^■«'*'*v**"v ^Vft '*^ r ' ^ ' V f 




Wavelength (E) 




Wavelength (A) 



Fig. B.6. Best-fit model for MPG 113 (in red) compared with observed spectra (in black): the STIS spectrum (upper left) and the FLAMES+ 
LIVES spectra (upper right). The bottom plot shows a comparison with the LIVES spectrum for spectral regions around O in /15592 (left), C m 
^5696 (middle), and He i ,15876 (right). 



photospheric absorption, and the upper limit on the mass-loss without degrading the fit to C in /U176. The weak C in AA464-7- 
rate is M — 6x 10~ 9 M yr~', 50-52 supports the low value for the carbon abundance. 



Depending on the v rot used for stellar evolution models, the 
location of MPG 523 in the H-R diagram indicates contradictory 
ages, either very young (with an age younger than 1 Myr) or 
much older (with an age of about 5 Myr). On the other hand, the 
carbon and nitrogen surface abundances are characteristic of sig- 
nificantly more evolved status (carbon depleted and nitrogen en- 
riched). The N/C ratio is one of the highest of the whole sample. 
We note that the C iv Al 169 is too weak in our model, however, 
but could not be matched by increasing the carbon abundance 



Appendix B. 12: NGC 346-046 - 07 Vn 

The luminosity and effective temperature we derived for NGC 
346-046 indicates that this star is very close to the ZAMS, possi- 
bly leftward of the ZAMS within the error bars. Consistent with 
its spectral classification and with the results by Mokiem et al. 
(2006), this star has a very high rotation, the second highest of 
the sample, with usin i = 300 kms -1 . This strongly suggests 
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Fig. B.7. Left: Best-fit model for AzV 243 (red line) compared with the COS spectrum (black line). Right: Best-fit model for AzV 243 (red line) 
compared with the FLAMES spectrum (black line). 



that this evolutionary status of this star is to a large extent influ- 
enced by its rotation, possibly to the point of following chemi- 
cally homogeneous evolution. However, our modeling does not 
indicate any departure from the initial helium abundance, which 
is inconsistent with such an evolution (and with Mokiem et al. 
2006). The nitrogen-to-carbon ratio, though, is the highest of the 
whole sample. The C iv resonance doublet is present as a pure 
absorption profile, from which we derive an upper limit of M — 
6x 10~ 9 M yr~' for the mass-loss rate. 

Appendix B.13: NGC 346-031 - 08 Vz 

NGC 346-031 is the last member of the Vz class of our sample. 
We achieved a very good fit of both the UV and optical spectra 
and derived fundamental parameters such that its age is above 
3 Myr, in contradiction with the expected youth of Vz class- 
members. In contrast to Mokiem et al. (2006), we did not find 
it necessary to increase the helium abundance, but instead found 
that the very low mass-loss rate we derived can explain the spec- 
tral classification in this parameter range (see Sect. 5.2.2). 



Appendix B. 14: AzV 461 - 08 V 

The fundamental and wind parameters we derived for AzV 461 
are close to those of NGC 346-031, but for the micro-turbulence 
and rotational velocity. For the effective temperature we derived 
from the UV analysis (see Sect. 3.2), He i ,14471 and He n ,14686 
are too weak, while He n ,14200 and He n /14542 are well repro- 
duced. Since the optical (2dF) spectrum is at lower resolution 
and S/N, we mostly relied on the COS spectrum to derive the 
fundamental parameters. The very close T e ff and luminosities 
of AzV 461 and NGC 346-031 imply very close evolutionary 



status. Consistent with the higher osin i, we derived a higher 
nitrogen enrichment and carbon depletion in AzV 461. 

Appendix B.15: MPG 299 - 08 Vn / 09 V 

With ysin i - 360 kms~', MPG 299 is the fastest rotator of our 
sample. The analysis of the FLAMES spectrum yields parame- 
ters in good agreement within the error bars with those by Mok- 
iem et al. (2006). We confirm their conclusion that helium is 
enriched in MPG 299. Despite the significant centrifugal cor- 
rection to the surface gravity, we found that the spectroscopic 
mass of MPG 299 agrees very well with its evolutionary mass. 
Although the location of MPG 299 in the H-R diagram indicates 
that it is a young object, we note that evolutionary tracks with 
initial u ro , higher than 360 kms -1 should be considered in this 
specific case (see Sect. 5.2.1). For MPG 299, the nitrogen abun- 
dance is such that higher M would produce too strong a PCygni 
profile at the N v resonance line, which is not observed. 

Appendix B.16: MPG 487 - 08 V 

MPG 487 was originally studied in the FUV and optical in 
Bouret et al. (2003), then Heap et al. (2006). Subsequently, 
Massey et al. (2009) used their own optical spectrum to pro- 
vide a different analysis of this star. They argued that the cooler 
J e ff we derived (along with Heap et al. 2006) was caused by a 
wrong correction of the moonlight continuum contribution that 
contaminates the AAT spectrum (same comments hold for MPG 
368). Such a contribution was indeed present and mis-corrected 
in our initial study. On the other hand, the claim by Massey et al. 
(2009) that we were unable to obtain a satisfactory fit to the op- 
tical without assuming an unlikely low value for the metallicity, 
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Fig. B.8. Left: Best-fit model for AzV 446 (red line) compared with the COS spectrum (black line). 



is wrong, as we demonstrate here. The metallicity was derived 
on the basis of the STIS spectrum analysis (see also Bouret et al. 
2003), which is not affected by moonlight contamination. The 
derived value allows the best fit to the very weak lines of iron 
ions present on this 5775' spectrum. Assuming higher metallicity 
would lead to an even lower T e s than the one we derived, hence 



increasing the discrepancy with the Massey et al. (2009) results. 
Moreover, the value quoted by these authors for the projected 
rotation rate of MPG 487 (usin i - 160 kms~') simply cannot 
be true (a typo, we assume). The lines of this star, whether in the 
FUV or in the optical are really narrow, indicating a v sin i that 
must be 20 km s _1 or even lower (see e.g. the C m A\ 176 lines). 



Article number, page 31 of 43 



X 

E 

0J 

"o 

< 




4x10" 
2x10' 






Wavelength (E) 



E 

U 

O 

2 






3650 


3690 




' 


*f* 


j^iwU^^i. 




0.90 




1 - 




D.77 




. 


-a 
a 

N 


0.66 




iii 



3730 
1.02 



fV - vu 1 



5010 5020 5030 5040 5050 



5587 5592 5597 




J5 


5405 5415 


54:^ 


1 
1 


1 1 II 
1 1 II 



491049SO 6460 6500 

Wavelength (E) 




5860 5885 



Fig. B.9. Left: Best-fit model for MPG 356 (red line) compared with the STIS spectrum (black line). Right: Best-fit model for MPG 356 (red 
line) compared with the UVES spectrum (black line). 



The SED of this star furthermore suggests that MPG 487 is 
a binary with a rather late-type companion (cf. Sect. 5.1). This 
conclusion was already mentioned in Bouret et al. (2003), al- 
though it was not as firmly established as it is now. 

Finally, our new UVES spectrum of MPG 487 (cf. Sect. 2) 
does not suffer from moonlight contamination and has an ex- 
cellent S/N ratio. We still find the same value of T e ff - 35 kK 
that we derived in Bouret et al. (2003). The surface gravity is 
higher than in Bouret et al. (2003) and more compatible with the 
value from Massey et al. (2009); we used the full line-set from 
the Balmer series to constrain this parameter. We also derived a 
non-standard abundance pattern for this star, with nitrogen being 
depleted while carbon seems enhanced compared with the stan- 
dard chemical composition of HII regions of the SMC by Kurt 
&Dufour(1998). 



Appendix B.18: AzV 468 - Q8.5 V 



Fundamental parameters were derived for the COS spectrum, ex- 
cept for the surface gravity, which we adopted as log g - 4.0, 
because no reliable diagnostics are present in the UV range. Be- 
cause of the moderate usin i of this star, photospheric iron and 
carbon lines are clearly defined and were used to constrain r e jf 
to a good accuracy. The COS spectrum of AzV 468 shows no 
trace of wind, either at C iv or at N v resonance doublets. The 
derived mass-loss rate is indeed very low and is merely an upper 
limit. 



Appendix B.19: MPG 682 - 09 V 



Appendix B. 1 7: AzV 267 - 08 V 

With usin i - 220 kms -1 , AzV 267 is another fast rotator of 
the sample. Because we only had the UV COS spectrum for 
this star, we constrained the effective temperature from ratios 
between iron ions and C iv A 1169 to Cm A 1176. Although it 
is a later-type O star, the C iv UV resonance doublet presents a 
clear absorption extending to the blue that we were able to use to 
constrain the mass-loss rate and the clumping filling factor. The 
fit to the N v resonance doublet is very good for the derived wind 
parameters. The terminal velocity we adopted to fit the lines is 
lower than for the other stars of the same type (i.e., with no wind 
detection). 



The fit to the COS UVES and FLAMES spectra of MPG 682 is 
excellent, including the weak lines of N in AA4634--42. Our T^ 
is cooler by 2000 K than the value reported in Mokiem et al. 
(2006), while log L/L Q are very similar. The wings of the hydro- 
gen lines (including the Balmer series below 3860 A) indicate 
that \ogg = 4.1, which yields a spectroscopic mass close to 27 
M , significantly higher than the evolutionary mass we derived 
from the H-R diagram (« 21 M ). The surface gravity quoted 
in Mokiem et al. (2006) is even higher (4.2 dex), but consis- 
tent with ours, within the error bars. The absence of UV wind 
profiles confirm that this star is in the weak-wind regime, as ex- 
pected from its spectral type and luminosity class. Our mass-loss 
rate is about two orders of magnitude lower than the mass-loss 
rate in Mokiem et al. (2006). 
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Fig. B.10. Best-fit model for AzV 429 (red line) compared with the COS spectrum (black line). 



Appendix B.20: AzV 326 - 09 V line from UV to the lower resolution 2dF spectrum in the optical 

is very good and we are confident that the derived parameters 
AzV 326 is very likely a multiple star (see the discussion in Sect. provide a me aningful representation of the primary component 
5.1). However, even when we corrected for a possible binary of this (vet hypothetical) system 
component, there were some remaining disagreement in the fit 
to the COS SED and we finally chose to normalize the spectrum, 
which we present in Fig. B.20). The fit to the individual spectral 
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Fig. B.ll. Left: Best-fit model for MPG 523 (red line) compared with the STIS spectrum (black line). Right: Best-fit model for MPG 523 (red 
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Appendix B.21: AzV 189 - 09 V 

Here again, we are dealing with a potential binary. The correc- 
tion applied assuming that this system consists of two identical 
components yields an excellent fit to the whole SED from UV to 
NIR with a single parameter set. For the same parameters, the fit 
to the optical 2dF spectrum is also very good. This spectrum was 
used to constrain the surface gravity. For the derived log g the 
spectroscopic mass is higher than the evolutionary mass. This 
could indicate that the adopted luminosity is too low, although a 
change by 0. ldex in loggr would also reconcile the two masses 



and is not excluded given the lower S/N and spectral resolution 
of the 2dF spectrum. 



Appendix B.22: AzV 148 - 08.5 V 

AzV 148 is the third star with photometry and UV flux unusual 
for its spectral type, and which would correspond to a luminosity 
typical of 06 V The lower v sin i for this star allowed us to derive 
accurate T s $ from ionization ratios, despite the lack of optical 
spectrum. The value we derived is consistent with those of the 
stars with spectral type 09 V. Surface gravity, on the other hand, 
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Fig. B.15. Left: Best-fit model for MPG 299 (red line) compared with the COS spectrum (black line). Right: Best-fit model for MPG 299 (red 
line) compared with the FLAMES spectrum (black line). 



Article number, page 36 of 43 



Bouret et al.: Massive stars at low metallicity 




a 

E 

0) 

3 

"o 

< 



10"' 3 


1 


1 


1 


J 


,Q-« 
















,o- 


\ 














- 


10 '* 


\ 














- 





1 


J 


1 


1 


1 


1 





V if 

1 , , , , 1 


T: 



1200 


1242 


1283 


32S 


1367 


1408 14 






immml^U. 


1 


1 


1 


1 


10 ia 






IO Ha 














10"' 4 












" 









1 


1 


II 1 


1 




mmmfiy^^ 



5412 5416 



■y ^ ^( Sm t 1 * ' ^ 




1 ' ' 




6510 5530 6550 6570 6590 6610 6630 



Wavelength (A) 



Wavelength (A) 



Fig. B.16. Left: Best-fit model for MPG 487 (red line) compared with the COS spectrum (black line). Right: Best-fit model for MPG 487 (red 
line) compared with the UVES spectrum (black line). 

was adopted as log g - 4.0, which translates into a spectroscopic 
mass significantly higher than the evolutionary mass. Like in 
AzV 189, though, a change by 0. ldex in log g would be enough 
to reconcile the two estimates. The C iv /1/11548-1550 line shows 
some blue extension, blueward of the interstellar components. 
The profile, together with the N v resonance profiles, is well re- 
produced by adopting a mass-loss rate three times higher than 
for the other 09 V stars of the sample. 

Appendix B.23: MPG 012 - 09.5 V / B0 IV 

This star was previously studied in the UV and optical in Bouret 
et al. (2003) and in the optical in Mokiem et al. (2006). The 
higher S/N ratios of the UVES and FLAMES spectrum led us to 
revise log g (0.05 dex) slightly upward compared to Bouret et al. 
(2003), while T e g was unchanged. The fit to the three spectra 
is excellent except for He i /i5876; the synthetic profile is not as 
strong and broad as observed. This behavior cannot be explained 
at this point. The star is still the coolest and lowest gravity ob- 
ject of the whole sample. The fundamental parameters suggest a 
clearly evolved status for this object, as indicated by its age, the 
highest by far among the stars in NGC 346. The surface abun- 
dances are also typical of more advanced evolution, as noted in 
Bouret et al. (2003). Because it is at the outskirts of this H II 
region, this object might not be coeval with this region, as pre- 
viously suggested by Walborn & Howarth (2000), Bouret et al. 
(2003), and Mokiem et al. (2006). 
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Fig. B.17. Left: Best-fit model for AzV 267 (red line) compared with the COS spectrum (black line). 
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Fig. B.20. Left: Best-fit model for AzV 326 (red line) compared with the COS spectrum (black line). Right: Best-fit model for AzV 326 (red 
line) compared with the 2dF spectrum (black line). 
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Fig. B.21. Left: Best-fit model for AzV 189 (red line) compared with the COS spectrum (black line). Right: Best-fit model for AzV 189 (red 
line) compared with the 2dF spectrum (black line). 
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Fig. B.22. Left: Best-fit model for AzV 148 (red line) compared with the COS spectrum (black line). 
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Fig. B.23. Best-fit model for MPG 012 (in red) compared with observed spectra (in black): the STIS spectrum (upper left) and the FLAMES 
spectrum (upper right). The bottom plot shows a comparison with the UVES spectrum. 
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